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Abstract 

We perform a complete analysis of isospin breaking in K — > 2tt amplitudes in 
chiral perturbation theory, including both strong isospin violation {m u ^ rrid) and 
electromagnetic corrections to next-to-leading order in the low-energy expansion. 
The unknown chiral couplings are estimated at leading order in the 1/N C expansion. 
We study the impact of isospin breaking on CP conserving amplitudes and rescatter- 
ing phases. In particular, we extract the effective couplings gs and 527 from a fit to 
K — * tttt branching ratios, finding small deviations from the isospin-limit case. The 
ratio Re^4o/R e ^2 measuring the AI = 1/2 enhancement is found to decrease from 
22.2±0.1 in the isospin limit to 20.3±0.5 in the presence of isospin breaking. We also 
analyse the effect of isospin violation on the CP violation parameter e', finding a de- 
structive interference between three different sources of isospin violation. Within the 
uncertainties of large- N c estimates for the low-energy constants, the isospin violating 
correction for e' is below 15 %. 
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1 Introduction 

A systematic treatment of isospin violation in nonleptonic weak interactions is needed for 
many phenomenological applications. The generically small effects induced by electromag- 
netic corrections and by the quark mass difference m u — are enhanced in subdominant 
amplitudes with A J > 1/2 because of the AI = 1/2 rule. For one, a quantitative under- 
standing of the AI = 1/2 rule itself is only possible with isospin violating effects included. 
Another area of application is CP violation in the K° — K° system where isospin breaking 
is crucial for a precision calculation of e'/e. 

Isospin violation in K — > 2ir decays has already been addressed in recent works Pfl El 
El IH El El El IE] • In this paper, we reanalyse the K — > -nit decay amplitudes to perform 
a comprehensive study of all isospin violating effects to next-to-leading order in the low- 
energy expansion of the standard model. More precisely, we shall work to first order in a 
and in m u — md throughout, but to next-to- leading order in the chiral expansion. In view of 
the observed octet dominance of the nonleptonic weak interactions, we therefore calculate 
to 0(Ggp 4 , G 8 (m u — m d )p 2 ,e 2 G 8 p 2 ) and to (^(G^tP 4 ) for octet and 27-plet amplitudes, 
respectively. 

At this order, many a priori unknown low-energy constants (LECs) appear. With few 
exceptions to be discussed in Sec. El we adopt leading large- N c estimates for the LECs. 
The advantage is that we employ a systematic approximation scheme with solid theoretical 
foundation that can in principle be carried through beyond leading order. On the other 
hand, the importance of subleading large- N c effects is at present not known in general. 
We shall estimate the uncertainties by varying the two scales entering those estimates: 
the renormalization scale for evaluating Wilson coefficients (short-distance scale) and the 
chiral scale (long-distance scale) at which the large- N c results are supposed to apply. 

In performing electromagnetic corrections, a careful analysis of radiative events is nec- 
essary as emphasized in Ref. We shall perform such an analysis for the new KLOE 
measurement \9\ of the ratio T(Kg — > ir + TT~['-f})/r(Ks — > tt tt°) with a fully inclusive 
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7r + 7r~[7] final state. The KLOE result influences the phase difference xo — Xi °f the two 
isospin amplitudes strongly. Together with this phase difference, the effective weak octet 
and 27-plet couplings G$, G27 will be the primary output of our analysis. With that out- 
put, several phenomeno logical issues can be addressed such as the relation of the phases 
Xo, Xi t° the s-wave pion-pion scattering phase shifts or the impact of isospin breaking on 
e'/e. 

The content of the paper is as follows. In the subsequent section, we introduce the decay 
amplitudes and the relevant effective chiral Lagrangians. The amplitudes at leading order 
in the low-energy expansion are presented in Sec. El The amplitudes at next-to-leading 
order are investigated in Sec. |H distinguishing between tt° — 77 mixing and all other contri- 
butions arising at that order. The amplitudes are divided into various parts depending on 
the source of isospin violation. The local amplitudes of next-to-leading order are explicitly 
given here. Sec. El analyses the LECs at leading order in 1/N C . To determine weak and 
electroweak LECs for N c — > 00, one needs input for the strong [up to 0(p 6 )] and electro- 
magnetic couplings [up to 0(e 2 p 2 )], in addition to the relevant Wilson coefficients. We 
discuss to which extent the necessary information is available. The numerical calculations 
of the various amplitudes to next-to-leading order in chiral perturbation theory (CHPT) 
are presented in Sec. El Dispersive and absorptive components of the loop amplitudes are 
given together with CP-even and CP-odd parts of the local amplitudes. Those amplitudes 
are then used in Sec. 13 to extract the lowest-order nonleptonic couplings G 8 , G27 and the 
phase difference Xo ~ X2 from K — > tctc data. We compare those quantities at lowest and 
next-to-leading order, the latter with and without isospin violation included. With this 
information, we then analyse the relation of the phase difference xo ~ Xi to the correspond- 
ing difference of tttt phase shifts. In Sec. El we discuss isospin violating contributions to the 
parameter e' of direct CP violation in K° — > titt decays. Sec. El contains our conclusions. 
Various technical aspects are treated in several appendices: next-to-leading-order effective 
chiral Lagrangians; explicit loop amplitudes; an alternative convention for LECs of lowest 
order; details for the analysis of the phase difference. 



2 Nonleptonic weak interactions in CHPT 

In this section, we define our notation for the K — ■> ixn amplitudes and we introduce the 
relevant effective chiral Lagrangians. 
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2.1 Decay amplitudes 



Using the isospin decomposition of two-pion final states, we write the K — > tttt ampli- 
tudes in the charge basis in terms of three amplitudes 1 Aai that are generated by the AI 
component of the electroweak effective Hamiltonian in the limit of isospin conservation: 

A- = Ai /2 + (A/ 2 + A/2 

Ao = Ax/ 2 - V2 (A 3/ 2 + A/2) (2.1) 



/2 — v ^ i ^3/2 -r SI5/2 

3 / 2 
A+o = - (A/2 - -A/2 

In the standard model, the AI = 5/2 piece is absent in the isospin limit, thus reducing the 
number of independent amplitudes to two. Each amplitude A n has a dispersive (VispA n ) 
and an absorptive (Abs A n ) component. In order to carry out phenomenological applica- 
tions and to keep the notation as close as possible to the standard analysis in the isospin 
limit, we write 

Ae' M = A/2 

A 2 e^ = A 3/2 + A 5/2 (2.2) 
A+e^ = ^3/2-2/3^5/2, 

where we explicitly separate out the phases xi- I n the limit of CP conservation, the 
amplitudes Ao, A 2 and A are real and positive. In the isospin limit, the Aj are the 
standard isospin amplitudes and the phases xi are identified with the s-wave nn scattering 
phase shifts 5i(y/s = M K ). 

For the phenomenological analysis (see Sees. [7| and |B1), we therefore adopt the following 
parametrization of K — > tttt amplitudes: 

A + - = A e ixo + A e iX2 
V2 

Ao = A e ixo -V2A 2 e tX2 (2.3) 
Ao = ~A + e^ ■ 

This parametrization holds for the infrared-finite amplitudes where the Coulomb and in- 
frared parts (defined in Sec. EJ) have been removed from A+-. 

1 We shall use the invariant amplitudes A n defined as follows: 

i I dx£(x) s 



(1 I clx L* [x ) \ 
e J J \K) = i {2ir) A (Pf - Pi) X (-» A n ) 
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In the absence of electromagnetic interactions = and therefore A 2 = A£- To set 
the stage, we extract the isospin amplitudes A , A 2 and the phase difference Xo ~ X2 from 
a fit to the three K — > titt branching ratios [2 El] : 

A = (2.715 ± 0.005) ■ 10" 7 GeV 

A 2 = (1.225 ± 0.004) ■ 10" 8 GeV (2.4) 
X0-X2 = (48.6 ±2.6)°. 

These values hold in the isospin limit except that the physical pion masses have been used 
for phase space. The substantial reduction in the phase difference Xo ~ X2 (from about 58° 
during the past 25 years jTUj) is entirely due to the new KLOE measurement of the ratio 

t(k s -> vr+7r-(7))/r( J fs:5 -> ttV) p. 



2.2 Effective chiral Lagrangians 

In the presence of isospin violation, the physics of K — > tttt decays involves an interplay 
of the nonleptonic weak, the strong and the electromagnetic interactions. Consequently, 
a number of effective Lagrangians are needed to describe those transitions. We use the 
well-known Lagrangian for strong interactions to 0(p 6 ) [TT ] IT2* j IT3]. the nonleptonic weak 
Lagrangian to 0(GfP^) O EE1 EH1 E! , ^ ne electromagnetic Lagrangian to 0(e 2 p 2 ) [TBI HHj 
and, finally, the electroweak Lagrangian to 0(e 2 G 8 p 2 ) [2U1 12T1 122] . 

Only the leading-order (LO) Lagrangians are written down explicitly here. The relevant 
parts of the next-to-leading-order (NLO) Lagrangians can be found in App. |X] along with 
further details. 

Strong Lagrangian: 

Strong = ^-(D^UD^ + X U ] + X ] U) 

+ ^LiOf + J^XiF-'of. (2.5) 

i i 

F is the pion decay constant in the chiral limit, the SU(3) matrix field U contains the 
pseudoscalar fields and the scalar field x accounts for explicit chiral symmetry breaking 
through the quark masses m u ,md,m s . The relevant operators Of are listed in App. El 
The LECs Xj of 0(p 6 ) will only enter through the large- N c estimates of the electroweak 
couplings in Sec. El (A) denotes the SU(3) flavour trace of A. 

Nonleptonic weak Lagrangian: 

£ W eak = G 8 F 4 (XDWD^U) + G 27 F 4 (L^L^ + -L^L^ 

+ G s N t F 2 0^ + J2 G 27 D t F 2 27 + h.c. (2.6) 

i i 
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The matrix = ill^D^U represents the octet of V — A currents to lowest order in 
derivatives; A = (A 6 — 2A7) /2 projects onto the s — > d transition. Instead of G$, G27 we 
will also use the dimensionless couplings g$, g?j defined as 

Gf 

G 8 ,27 = V ud V us 98,27 ■ (2.7) 

One of the main tasks of this investigation will be the determination of gs, g27 in the 
presence of isospin violation to NLO. The LECs iVj, Di of 0(GfP 4 ) are dimensionless. The 
monomials Of, Of 7 relevant for K — > 2n transitions can be found in App. 1X1 

Electromagnetic Lagrangian: 

£ clm = e 2 ZF 4 (QU^QU) 

+ e 2 Y, K, t F 2 0f p2 . (2.8) 

i 

The quark charge matrix is given by Q = diag(2/3, —1/3,-1/3). The lowest-order LEC 
can be determined from the pion mass difference to be Z ~ 0.8. The NLO LECs Ki are 

2 2 

dimensionless and the relevant monomials 0\ p can again be found in App. lAl 
Electroweak Lagrangian: 

Cew = e 2 G 8 g ewk F G (XU^QU) 

+ e 2 J2 G 8 ZiF 4 0f w + h.c. (2.9) 

i 

The value of the LO coupling g ewk is discussed in Sec. El The LECs Z { are dimensionless 
and the associated monomials Of w are collected in App. El We do not include isospin 
violating corrections for 27-plet amplitudes. 

The low-energy couplings Li,Ni, Di, Ki, Zi are in general divergent. They absorb the 
divergences appearing in the one-loop graphs via the renormalization 



(2.10) 



where v x is the chiral renormalization scale and the divergence is included in the factor 



Li 


= L\{ 




+ r i A(v x ) 




= Nl 






A 


= DI 




+ diA(v x ] 


Ki 


= K 'i 




+ Ki k(v x 


Zi 


= zk 


v x) 


+ ZiA(v x ) 



r 1 1 



log(4vr)+r / (l) + l 



:2.1V, 



XJ (4vr) 2 Id- 4 2 

The divergent parts of the couplings are all known ^TJ Q21 EH 122] and they allow for 
a nontrivial check of the loop calculation. On the other hand, many of the renormalized 
LECs contributing to the decay amplitudes are not known. Our strategy will be to use LO 
large- N c estimates. A comprehensive discussion of all relevant LECs will be presented in 
Sec. El 
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3 Amplitudes at leading order [O(GFp 2 1 e 2 Ggp )] 

With the most general effective chiral Lagrangian of the previous section, we can now 
proceed with the construction of physical amplitudes. At LO in the low-energy expansion, 
the procedure is straightforward: chiral power counting tells us that the amplitudes are 
obtained by summing all tree-level Feynman diagrams with one insertion from either £ wea k 
of 0{GfV 2 ) or £ew of O(e 2 G$p ), at most one insertion of £ e im of O(e 2 p ) and any number 
of insertions from the 0(p 2 ) part of the strong Lagrangian (|2.5J1 . 

In addition to contributions proportional to the electroweak coupling g e wk, isospin 
breaking occurs also in the pseudoscalar mass matrix, generating in particular non-diagonal 
terms in the fields (^3, 7Tg) (tt° — 77 mixing). Upon diagonalizing the tree-level mass matrix 
one obtains the relation between the LO mass eigenfields (tt°, 77) and the original fields 
(7r 3 , 7T 8 ) (to first order in m u — rrid): 



( 7T 3 \ ( 1 -e® 



( 



5 



(2) I 



7T° 



V v 



(3.1) 

LO 



with the tree-level tt° — 7] mixing angle given by 



£ ( 2 ) = V3 m d -m u ^ 
4 m s — m ' 

where m stands for the mean value of the light quark masses, 

rh = ^(m u + m d ) ■ (3.3) 

The physical amplitudes are then obtained by considering the relevant Feynman graphs 
with insertions from the LO effective Lagrangian expressed in terms of the LO mass eigen- 
fields. 

Apart from 7r° — 77 mixing, isospin breaking manifests itself also in the mass differences 
between charged and neutral mesons, due to both the light quark mass difference and 
electromagnetic contributions. We choose to express all masses in terms of those of the 
neutral kaon and pion (denoted from now on as Mx and M^, respectively). In terms of 
quark masses and LO couplings (B is related to the quark condensate in the chiral limit 
by (0|gg|0) = —F 2 B ), the pseudoscalar meson masses read: 



(3.4) 

, - to) + 2 e 2 ZF 2 
—= B (m s - m) . 



Ml 


= 2B m 




= Ml + 2e 2 ZF 2 


Ml 


= B (m s + m d ) 


M 2 K± 


= Ml - ^_ Bo( 


Ml 


= \UMl-Ml) 



7 



We are now in the position to write down the three independent amplitudes relevant 
for K — > mi decays. In the physical "charge" basis the LO amplitudes are 



A+- 
-4-00 
^4+0 



-V2G 27 F (m 2 k 



Ml 



V2G 8 F 



M l K -Mi- e l F l (g ewk + 2Z) 



-V2G 27 F (M 2 K - Ml) + V2G 8 F (M 2 K - Ml) f 1 " 4=< 



■(2) 



(3.5) 



~ G 27 F ( M 2 K 



Ml ) + G$F 



(Ml - Ml) A - e 2 F 2 (g cwk + 2Z) 



We recall that we do not include isospin violation for the 27-plet amplitudes. In the isospin 
basis, more convenient for phenomenological applications, the LO amplitudes are given by 
(see Eq. (|2.1j) for the relation between the two bases) 



A 



1/2 



^G 27 F(M 2 K -Ml) 
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+ V2G 8 F 



(M 



M: 



A 3 / 2 
A 5 / 2 



10 
. 



(3.6) 



G 27 F (Ml - Ml) + G S F (Ml - M$) ^= - ± e 2 F 2 (g cwk + 2Z) 



The parameter F can be identified with the pion decay constant F n at this order. The 
effect of strong isospin breaking (proportional to e^) is entirely due to 7r° — 7] mixing at 
LO. Electromagnetic interactions contribute through mass splitting (terms proportional 
to Z) and insertions of g ew k- As a consequence of imposing CPS symmetry (22 on the 
effective Lagrangian, electromagnetic corrections to the octet weak Hamiltonian do not 
generate a AI = 5/2 amplitude at LO in the quark mass expansion. 



4 Amplitudes at next-to- leading order [0(GfP A , e 2 Ggp 2 )] 

Let us now outline the construction of NLO amplitudes. As always, chiral power counting is 
the guiding principle: it tells us that both one-loop and tree-level diagrams now contribute. 
In the one-loop diagrams, one has to consider one insertion from either £ wea k of 0{GfP 2 ) or 
£ew of O(e 2 Gsp ), at most one insertion of £ e im of O(e 2 p ) and any number of insertions 
from the 0(p 2 ) part of the strong Lagrangian, with the LO effective Lagrangians expressed 
in terms of the LO mass eigenfields. For the tree-level diagrams, one has to apply one 
insertion from the NLO effective Lagrangian and any number of insertions from the LO 
Lagrangian. The strangeness changing vertex can come from either the LO or NLO effective 
Lagrangians. This straightforward prescription leads to a large number of explicit diagrams 
for each mode, due to several topologies and several possibilities to insert isospin breaking 



S 



vertices from the LO effective Lagrangian (in the weak vertex, in the strong vertex, in 
the internal propagators, in the external legs). We begin with the well-defined class of 
NLO corrections to the pseudoscalar meson propagators, focusing afterwards on the other 
corrections. 



4.1 



7T 



77 mixing at NLO 



As in the LO case, it is convenient to first analyse isospin breaking in the two-point 
functions (inverse propagators) and to define renormalized fields in which the propagator 
has a diagonal form with unit residues at the poles 2 . At NLO two main new features arise: 

• Not only the (713, 7r§) mass matrix acquires off-diagonal matrix elements, but also the 
purely kinetic part of the propagator does so. 

• Electromagnetic interactions contribute to this phenomenon, in addition to the up- 
down mass splitting. 

Results on NLO mixing effects induced by quark mass splitting already appear in Refs. (TTJ 
I2HIE] 5 wn il e electromagnetically induced effects were considered in [23122] ■ Here we follow 
the formalism outlined in Ref. |24j . treating strong and electromagnetic effects simulta- 
neously. In the LO mass eigenfield basis, the NLO inverse propagator (an eight-by-eight 
matrix) can be written as follows: 



n(g 2 ) 



g 2 i 

eg 2 



M 2 - fl(g 2 ) 



(4.1) 



where M 2 is the diagonal LO mass matrix and C, D are symmetric matrices, which are 
diagonal except for their restriction to the (7r°, 77) subspace. The relation between the LO 
and NLO mass eigenfields (collected in a vector <p a ) can be summarized as follows: 




NLO 



(4.2) 



where W is an antisymmetric matrix, non-vanishing only in the (77°, 77) subspace. Except 
for the reduction to this subspace, Eq. (J4.2)) is just the familiar field renormalization, with 
wave function renormalization given by Zi = 1 + Cu. Focusing on the (it , rf) sector, we note 
that W is characterized by a single entry called (23- This quantity is UV finite and 
represents a natural generalization to 0(p 4 ) of the tree-level mixing angle e^ 2 \ Explicitly, 
the relation between the (tt°, 77) mass eigenfields at LO and NLO is given by 



\ 1 



( ^ C w n 



LO 



-e 



r \ 

(4) , 



+ 



£ (4) 




(4.3) 



NLO 



2 This way, no further wave function renormalization effect has to be included. 
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Eqs. (13.1(1 and (|4.3(1 give the full relation between the original fields (7^, 7rg) and the NLO 
mass eigenfields. We do not report here the factors (7 a &, as they are UV divergent and 
make sense only in combination with other terms in the full amplitudes. We do report, 
however, the expression for because the replacement —> e^ 2 ' + in Eqs. (|3. 5 J3.fi J) 
gives rise to a well defined (UV finite) subset of the NLO corrections. Breaking up 
into contributions from strong (S) and electromagnetic (EM) isospin breaking, one gets 



r(4) 



.(4) 



e (4) + e (4) 

fc S ^ fc EM 



2e P) 



3(4vrF) 2 (M2 - Ml) 



2\2 



(47r) 2 64[3L 7 + L^ x )] (M 2 K - M>) 



Ml 



Ml 



- M 2 {M 2 K - Ml) log -n- + M>(M 2 K - 3M 2 ) log 



-(4) 
-EM 



2M 2 K {M 2 K -2Ml)\og 



M 



K 



v- 



2M 2 K (M 2 K - Ml] 



2^3 



108?r(M2 -Ml) 



9M Z K Z log 



+ 1 



//f 



+ 2M 2 K (4n) 2 2U;(u x ) + 3U r 3 (^ 



+ mKatt) 2 2u;{v x ) + 3u;(u x ) - w. 



(4.4) 



The electromagnetic LECs U{ are linear combinations of the K\ (defined in Sec. 15.5(1 . 



4.2 Remaining NLO contributions: a guided tour 

Having dealt with the propagator corrections in the previous section, we now describe the 
remaining contributions to the K — > ttti amplitudes at NLO, starting with the one-loop 
terms. There are two main classes of contributions: loops involving only pseudoscalar 
mesons (Fig. HJ and loops involving virtual photons (Fig. Hj). In the isospin limit, contri- 
butions to the amplitudes arise from the topologies of Fig. ^ by inserting the LO weak 
vertices proportional to G 8 or G 27 in the Lagrangians (|2.fi() and ()2.9|) . Given the large 
suppression of G27/GS, we consider in this work only isospin breaking effects generated 
through the octet component of the effective Lagrangian. We are therefore interested 
in the terms proportional to e^Gs (strong isospin breaking) and e 2 Gs (electromagnetic 
isospin breaking). 

Strong isospin breaking terms (e^Gs) at NLO come from several sources: 

• Explicit terms ~ (m u — rrid) in the strong vertices of Fig. ^ obtained by expressing 
^strong in terms of the LO mass eigenfields. 

• Mass corrections in the internal propagators, for which we use the LO diagonal form 
(and the corresponding mass relations of Eq. (|3.4(l ). 
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Figure V. Topologies for purely mesonic loop diagrams contributing to K — > txtx: the filled 
circles indicate AS = 1 vertices of lowest order. Wave function renormalization diagrams 
are not shown. 



• Mass corrections arising when external momenta are taken on-shell (using again 
Eq. flUD)- 

The combination of these effects leads in principle to non-linear contributions in the isospin 
breaking parameter. We have chosen to expand the final expressions for the amplitudes to 
first order in . 

Electromagnetic isospin breaking terms (e 2 G$) at NLO can be naturally divided into 
three categories: 

• e 2 ZG$: these arise exactly in the same way as the strong isospin breaking terms (see 
discussion above). 

• e 2 G 8 g ew \ l : these arise from insertions of the g ew k vertices of £ew m the topologies 
of Fig. ^ keeping all other contributions (masses and strong vertices) in the isospin 
limit. 

• e 2 G§: these arise from the photonic diagrams of Fig. |21 using the LO weak vertices 
of £ wea k proportional to Gg. This class of contributions to A + - is infrared divergent. 
We regulate the infrared divergence by means of a fictitious photon mass M 7 . The 
cancellation of infrared divergences only happens when one considers an inclusive sum 
of K — > 7T7T and K — > infy decay rates and we postpone details on this point to Sec.[7| 
At this stage, we split the photonic correction to A + - into an "infrared component" 
A^_(Mry) (to be treated in combination with real photons) and a structure dependent 

part ^4+1, which is infrared finite and has to be used together with the non-photonic 
amplitudes in Eq. (j2.1|) . Clearly, an arbitrary choice appears here as one can shift 
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Figure 2: Topologies for meson-photon loop diagrams contributing to K — > txtx: the filled 
circles indicate AS = 1 octet vertices of lowest order. 



infrared finite terms from A+L to A l +_(M^). This also implies that the isospin 
amplitudes all depend on this choice. The observables, however, are only affected by 
this ambiguity at order a 2 . Azi(M 7 ) has the following structure: 



in terms of the function 5 + _(M 7 ) reported in App. El 

This concludes our description of one-loop contributions to K — ► ttti amplitudes. The 
NLO local contributions arise from tree-level graphs with insertions of one NLO vertex and 
any number of LO vertices, according to the topologies depicted in Fig. El 



Figure 3: Diagrams for NLO local contributions: the filled square denotes a NLO vertex. 



4.3 Structure of the amplitudes 

Having identified the various diagrammatic contributions to the physical amplitudes, we 
now introduce a general parametrization that explicitly separates isospin conserving and 
isospin breaking parts and allows to keep track of the various sources of isospin break- 
ing. Let n be the label for any amplitude. Then, including the leading isospin breaking 



Af_ (M 7 ) = V2G 8 F (M 2 K - Ml) a B+_(M 7 ) , 



(4.5) 
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corrections (proportional to Gg), one has: 



A n = g 27 fJm 2 k -m 2 \a (27) 



(4.6) 



+ g s fJ (m 2 k -m 2 



A^ + e& Me) 



-e 2 F 2 



A^ + Z A n ^ + 5Wk A^ 



The meaning of the amplitudes A^' can be inferred from the superscript X. A^\A^ 7 ' 
represent the octet and 27-plet amplitudes in the isospin limit. Affl represents the effect of 
strong isospin breaking, while the electromagnetic contribution is split into a part induced 
by photon loops A^ and the parts induced by insertions of Z and g ew k vertices {A^' and 
A%\ respectively). 

At the order we are working, each of the amplitudes A n xs> can be decomposed as follows: 



A (x) 



1 + A L A^ + A c A n x \ 
A L A™ + &c AW 



if ) ^ 



(4.7) 



if 



il x) = 



with 



A L A^ 
&cA^ 



LO contribution [Eqs. (j3.5|3.6j) 
NLO loop correction 
NLO local correction . 



The amplitudes a n x \ A L A n x ^ and A c A n x ^ are dimensionless and we have chosen to 
normalize the NLO contributions to the LO contributions whenever possible. Moreover, in 
Eq. (|4.6|) we have traded the constant F for F n , the physical pion decay constant at NLO. 
The relation between the two is given explicitly by flj |2Zj 



1 - 



F 2 



L\(y x ) (Ml + 2M 2 K )+Ll{v x )Ml 



+ 



2(4tt) 2 F 2 
=V (Ml -Ml 

Vs y K 



M 2 M 2 
2M>g^ + M|log M * 



8Ll(u x ) 



vt 



F 2 2(4vr) 2 F 2 



1 + log 



M£ 



(4.8) 



Both A^A n x ^ and AqA^ individually are UV divergent and scale dependent. Only in 
their sum the UV divergence and the scale dependence cancel, providing a valuable check 
on the calculation. The explicit form of the various loop contributions is given in App. [B] 
while the local amplitudes are reported in the next subsection. 



13 



4.4 Local amplitudes 



The NLO K — > titt local amplitudes receive contributions from the NLO couplings Lj, iVj, 
Dj, i^j, Zi in the effective Lagrangians of Sec. 121 Following Ref. [2H], it is convenient to 
define the combinations 

Ml M 2 
A c = --^(4^ + 32^) --L(ULl + 16Ll) 

Ag w) = -^(4L5 + 48L;)-^(20L5 + 24LJ) . (4.9) 

In terms of the quantities defined above, the finite parts of the NLO local amplitudes have 
the form reported below. In this section we use the notation D i: Ni, Zi as a shorthand for 
the ratios of NLO to LO chiral couplings [g B Dt)/g B , (g 8 Ni)/g s , (g 8 Zi)/g 8 . 



AI = 1/2 amplitudes 



A 



A A® 



A C A 



1/2 



Ac4] 2 



+ 



+ 
+ 



2M 2 

A c - 
2 Ml 

A c - 

2 M 2 K 

2M 2 



(D r A -D r 5 -9D r 6 + 4D r 7 ) 



(-6D[-2D r 2 + 2Dl + 6D r 6 + D r 7 ) 
2 M 2 

-yf (-N r 5 +2N r 7 -2N r 8 -N r 9 ) 
(-2 N£ - 4 N; - Nl + 2 N[ + 4 A+, +2N- 



<12J 



{M 2 K - Ml 
F 2 



(96LI + 32 (3Z/ 7 + Lg)) 



(iV 5 r + 6 Nq + 12 A/7 — 8 A/g — A/g - 4 A^[ - 8 AT[ 2 - 12 A/ 



+ — ^ (14 A^ 5 r + 6 A/"g + 24 A/7 — 5 N£ — 26 N[ - 24 N r u - 10 N{ 2 - 12 AT[ 3 ) 



F 2 

~ (ew) _4M£ 



T (ew) 

2y/2 



F 2 



(2JV?-JVJ-JVJ) + 



2M 2 

7T 



(2 A/g + 4 A/^ + A/"g 



F 2 



(M 2 K - Ml 



(6 U[ + 4 f/ 2 r + C/ 3 r 



-^(36^ + 22^ + 3^ + 2 f/1 ) 



6F 2 



ZI + 24 ZI - 9 Z£ - 6 Z\ + 3 Z s r + 3 ZI + 2 Z 



2Z 



11 



Mi 



Ml 



(2^ + 4^) + ^ + 



2 Z[ 2 ) 



(4.10) 
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AI = 3/2 amplitudes 



A 

A C vA 3/2 



+ 
+ 



A C A$ 2 




A c - 
2M 2 K 






F 2 




+ 


2 Ml 
F 2 


A c Al% 




Ag w) 




+ 


Ml 
5F 2 


A c Ai% 




A( ew ) 


A C A$ 2 




2 
3 



2 M 2 

(Dl -Dl + AD r 7 ) + ^(-2D r 2 + 2Dl + D r 7 ) 



F 2 

32 (M| - M 2 



F 2 



yf (n; + qNq -2n;- n; -an[ q -%n[ 2 - i2N[,) 

\M 2 

(2N; + 6 Nq + N r 8 -2N r 10 - 10N{ 2 - 12 N r 13 ) 



+ 



M 2 



K 



5F 2 



(12 iV 5 r - 16 Nj + 20 iV 8 r + 8 iV 9 r ) 



(8 N r 5 + 16 Nj + 10 N r 8 + 12 A/"g 



F 2 5 Us F 2 {3 U2 + 5 U3 3 U4 J 

24 8 8 

-±z; + -zi-3z;-3Z r 9 -2Z[ --z r 11 --z 

Ml M 2 
-jg (2Z[ + AZ r 2 -Z r 6 ) + T ^(AZ[ + 2Z r 2 ) 



{M 2 K - Ml 
3F 2 



(4.11) 



AI = 5/2 amplitudes 



A c Af> 

l(7) 



A c 4 



4 
3 
2 
3 



\M 2 K -Ml) 



5F 2 



\M\ 



Ml) 



15 F 2 



{-12 Nl - 2AN r 7 + 12 N r y) 



;-18 f/ 3 r + 36 Z\ + 18 Z[j + 18 Z[ 2 ) 



(4.12) 



5 LECs at leading order in 1/N C 

Owing to the presence of very different mass scales (M n < Mk < A x Mw), the gluonic 
corrections to the underlying flavour- changing transition are amplified by large logarithms. 
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The short-distance logarithmic corrections can be summed up with the use of the operator 
product expansion and the renormalization group (30] > all the way down to scales 
Z^sd < rn c . One gets in this way an effective AS = 1 Lagrangian, defined in the three- 
flavour theory |HU E2 ESI E3 , 



Ceff 1 = ~^ V u d C E C^SD) Q^SD) , (5.1) 



which is a sum of local four-fermion operators Qi, constructed with the light degrees of 
freedom (m < //sd), modulated by Wilson coefficients C^sd) which are functions of the 
heavy masses (M > Hsb) and CKM parameters: 

Ci(fM S T>) = ^(a*sd) + ryi(nsT>) (5.2) 
V td V* 



T 



The low-energy electroweak chiral Lagrangian arises from the bosonization of the short- 
distance Lagrangian (|5.1|) below the chiral symmetry breaking scale A x . Chiral symme- 
try fixes the allowed operators, at a given order in momenta, but the calculation of the 
corresponding CHPT couplings is a difficult non-perturbative dynamical question, which 
requires to perform the matching between the two effective field theories. 

The 1/N C expansion provides a systematic approximation scheme to this problem. At 
leading order in 1/N C the matching between the three-flavour quark theory and CHPT 
can be done exactly because the T-product of two colour- singlet quark currents factorizes. 
Since quark currents have well-known realizations in CHPT the hadronization of the weak 
operators Qi can then be done in a straightforward way. As a result, the electroweak chiral 
couplings depend upon strong and electromagnetic low-energy constants of order p 2 ,p 4 ,p 6 
and e 2 p 2 , respectively. 

5.1 Weak couplings of 0(G F p 2 ), O(e 2 G 8 p ) 

At lowest-order [0(GfP 2 ), O(e 2 Gsp )}, the chiral couplings of the nonleptonic electroweak 
Lagrangians (|2.6|) and ()2.9|) have the following large- N c values: 

2 3 

08° = -g Ci(a*sd) + - C 2 (/i S d) + C 4 (/xsd) - 16 L b B(/x S d) C 6 (// S d) 

9?7 = ^[Cl(/iSD)+C 2 ( yUS D)] (5.3) 

16 

(e 2 g 8 gewk)°° = -3 B(hsd) C s {i^sb) - y #(a*sd) C 6 (//sd) e 2 (K 9 - 2K W ) . 

The operators Qi (i ^ 6,8) factorize into products of left- and right-handed vector 
currents, which are renormalization-invariant quantities. Thus, the large- A^ factorization 
of these operators does not generate any scale dependence. The only anomalous dimensions 
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that survive for N c — > oo are the ones corresponding to Qq and Qs [13]. These operators 
factorize into colour- singlet scalar and pseudoscalar currents, which are /xsd dependent. 
The CHPT evaluation of the scalar and pseudoscalar currents provides, of course, the 
right /xsd dependence, since only physical observables can be realized in the low-energy 
theory. What one actually finds is the chiral realization of the renormalization-invariant 
products m q 9(1,75)?. This generates the factors 



SDj 




Ml 



1 2 



(m s + m d )(fj, S T)) 



1 



16M| 
F 2 



(2L 8 - Lj 



8M^ 

+ -^2-^5 



+ 



j(2M| + jg) 



(3L 4 - 4L e 



(5.4) 



in Eq. ()5.3jh which exactly cancel (S3 E3 EE E23 the /xsd dependence of 6*6,8 (a*sd) at 
large iV c . There remains a dependence at next-to-leading order. 

Explicitly, the large- N c expressions imply 3 

g? = (l.lO±0.05 (/ , SD) ±0.08 {L5) ±0.05 (ms) 
+ r f 0.55 d= 0.15 



f/'3D) ± - 20 (i 5 ) -ai6( ms ) 



921 



0.46±0.01 (AtSD) 



(5.5) 



{9s 9 



ewkj 



r 21.7 ±4.5 



(msd) 



± 1.0 



-9.1 



where the first uncertainty has been estimated by varying the renormalization scale /xsd 
between 0.77 and 1.3 GeV, the second one reflects the error on the strong LECs of order 
p 4 and e 2 p 2 , and the third indicates the uncertainty induced by m s jUJl which has been 
taken in the range [2H] {m s + m d )(fx SD = lGeV) = (156 ± 25) MeV. While the CP-odd 
component of g ew k is dominated by the electroweak penguin contribution (proportional to 
r l/8(/ i SD)), the CP-even part receives contributions of similar size from both strong (Q 6 ) 
and electroweak (Qs) penguin operators. Its large uncertainty within this approach reflects 
the GIM mechanism (^(/-^sd > m c) — 0). For the CP-even component, there exists an 
independent estimate, consistent with the one given here within the large uncertainties: 



Re (g 8 g 



cwk I 



-0.99 ±0.30 U 
-1.24 ± 0.77 (MSD) ± 0A0 (L5tKi) Eq. 



(5.6) 



In this work we shall always use the latter value, in order to perform a consistent analysis 
at leading order in 1/N C . 

3 According to the discussion presented in the following subsections, we use here L^(M p ) = (1.0 ± 0.3) • 
10-3 and [K% - 2K{ )(M p ) = -(9.3 ± 4.6) ■ lO" 3 . 
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Finally, the large-X c matching also produces the so-called weak mass term (see Ap- 
pendix El for notations): 

C wmt = G , 8 F 4 (\ x l l)+ h.c. , (5.7) 

with 

Gp 

(gX = -16 (L s + l -H 2 ^ B(fi SD )C 6 (tisD) ■ (5.8) 

We eliminate this term with an appropriate field redefinition [231 EH E3 , of the form 

U -> e ia Ue i(3 , (5.9) 

where the chiral rotation parameters (a and j3) are proportional to G' s . When applied 
to the strong effective Lagrangians of order p 4 and e 2 p 2 , the above redefinition generates 
monomials of the NLO Lagrangians of order G 8 p 4 and e 2 G$p 2 . The corresponding con- 
tributions to the couplings ggNi (of the form L n x (L 8 + 1/2 H 2 )) and g$Zi (of the form 
K n x (Lg + 1/2 H 2 )) need to be added to the results obtained by direct matching at large- 
X c . The complete results (reported in the next section) are independent of the unphysical 
LEC H 2 olO(p 4 ) HI]. 



5.2 Weak couplings of 0(G F p A ), 0(e 2 G 8 p 2 ) 

The large- N c matching at the next-to- leading chiral order fixes the couplings G$Ni, G 27 Di 
and GgZi of the nonleptonic weak and electroweak Lagrangians ()2.6|) and (|2.9|) . The oper- 
ators Q3 and Q5 start to contribute at 0(GfP 4 ), while the electroweak penguin operators 
Qi, Qg and Q10 make their first contributions at 0(e 2 Ggp 2 ). The contributions from the 
operator Q 6 at 0{G 8 p 4 ) involve the strong CHPT Lagrangian of 0(p 6 ) [13] (to avoid con- 
fusion with the Wilson coefficients Cj, the corresponding dimensionless couplings [32] are 
denoted here as Xj). With the definitions 

2 3 

Ci(/*sd) = --Ci(^sd) + gC 2 (AisD) + C 4 (//sd) 

3 2 

C 2 (/Usd) = +gC?i(/isD) - ^O^sd) + C 3 (jisd) ~ Csfjiso) , (5.10) 
the non-vanishing couplings contributing to K — > hit amplitudes are: 
(g 27 D 4 ) = 4L 5 g% 

(g 8 N 5 ) = -2L5C 1 ( / usD) + C 6 (/isD)5( /U sD)[-16X 1 4 + 32X 17 - 24X 38 - 4X 91 ] 
(g 8 N 6 ) = 4L 5 (5 1 (/i SD )+C 6 (// SD )S(// SD )[-32X 17 - 32X 18 + 32X37 + 16X 38 ] 
(g 8 N 7 ) = 2L 5 C 1 (/isD)+C , 6 (/iSD)5( / isD)[-32X 16 - 16X 17 + 8 X 38 ] 
(g 8 N 8 ) = AL 5 C 1 (fisD)+C 6 (fi S D)B(fi SD )[-16X 15 -32X 17 + lQX 38 } 

18 



(g 8 N 10 ) 

(9s N u ) 

(g 8 N 12 ) 
(g 8 N 13 ) 



C 6 (/x S d) B(hsv) [-64 L 5 L 8 - 8 X 34 + 8 X 38 + 4 X 91 ] 
C 6 (fi SB ) B(fi SD ) [-48 X 19 - 8 X 38 - 2 X 91 - 4 X 94 ] 
C 6 (// SD ) 5(// SD ) [-32X 20 + 4X 94 ] 
C 6 (^ S d) 5(^ S d) [128L 8 L 8 + 16X 12 - 16X 31 + 8X 38 



(5.11) 



2X 



4X, 



91J 



C 6 (/xsd) B(/i sb ) 



32 

256L 7 L 8 - y #12 



16 X 33 + 16 X 37 + -X 91 + 4 X 94 



Bosonization of the four-quark operators Qi in (|5.1|) leads to the following expressions 4 for 
the LECs Zt\ 



(9s Zx) 

(9s Z 2 ) 

(9s Z 3 ) 
(9s Z 4 ) 

(98 Z s ) 
(98 Z 6 ) 

(98 Z 7 ) 
(98 Z 8 ) 

(98 Zg) 

(9s Z w ) 

(9s Z n ) 
(9s Z 12 ) 

(98 Zxi) 

(98 Z\a) 
(9s Z lh ) 



SD; 



K a +64C 6 (/i SD )B(/is D )i! 



5 + 



K 



ii 



24 Cs(jisd)B(^sd) l 
e 2 8 

|c 1 (msd)^13-|64C 6 (msd)S( MS d) (#10 + #ll) ^8 

CiifMso) #13 - 64C 6 (/i SD ) S(^sd) (#10 + #ll) ^8 
-C^sd) #13 + 64C 6 (^sd) B(fi SB ) L 8 (X 10 + K n ) 

4 ~ , , „ „ „ „ „ x 64 
3 



d(MsD) (4X 1 + 3X 5 + 3X 12 ) 
|(X 5 + X 6 ) + 2(X 
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3 

#13 



C 6 (/i S D)5(/i SD ) (2K 7 + K 9 ) L 5 + 



Cio(/^sd) 



~ CeCAiso) S(a*sd) (#9 + #io + 3X n ) L 5 - 12 ^sd) ^sd) L 5 



(5.12) 



Ci(//sd) (8 #2 + 6 X 6 - 4X 13 ) - 32 C 6 (// SD ) B(// SD ) (2X 8 + X 10 
Ci(a*sd 

4 



X 



n 



4 ~ ~ 3 

o (Ci^sd) + C 2 (^sd)) #5 + — (O>0«sd) + C 10 (^sd)) 



2e 2 



(Ci(a*sd) + CaOW 

-2 CM/iSD) #13 + 4 (CM/lso) + C 2 (/i S D)) #6 



--Cj^SD) (#4 + #12 + #13) + ;:)> 



K R - 



2C 1 (jjl S d) (#4 + #13) 

-4C 1 (/i S D)#3 

2 ~ 64 

- Ci(psd) (#5 - #12 - #13) + y C b (/xsd) 5(^ S d) (#10 + #11) 

^(/iso) (-2X 6 + 4X 13 ) - 32C 6 (// S d) B(/x sd ) (X 10 + X n ) L 5 
dC/xsD) (-2X 6 + 4X 13 )-32C 6 (^ SD )S( ) u SD ) (X 10 + X n ) L 5 



3C 7 (A* 



SDj 



# 



We recall here that a matching ambiguity arises when working to next-to-leading order in 
the chiral expansion and at leading order in 1/X C : we cannot identify at which value of the 
chiral renormalization scale v x the large- X c estimates for the LECs apply. This turns out 

Z13, Z14, Z15 do not contribute to K — ► tttt amplitudes. 
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to be a major uncertainty in this approach. In order to account for this uncertainty, we 
vary the chiral renormalization scale between 0.6 and 1 GeV. The corresponding changes 
in the amplitudes are sub- leading effects in 1/N C and we take them as indication of the 
uncertainty associated with working at leading order in 1/N C . 

Finally, from the above expressions we see that in order to estimate the weak NLO 
LECs at leading order in 1/N C , one requires as input several combinations of strong LECs 
of order p 4 ,p e and e 2 p 2 . Below we summarize our knowledge of the needed parameters. 



5.3 Strong couplings of 0(p A ) 

It is well known that the limit N c —>■ oo provides an excellent description of the 0(p 4 ) 
CHPT couplings at v x ~ M p ^Hj- The leading-order contribution of Qq involves the 
LEC L5. The large- iV c value of this coupling can be estimated from resonance exchange 
|18j . Within the single-resonance approximation (SRA) Hlj . taking F = F n and 
M s = 1.48 GeV [45], one finds Lf = F 2 /(4M|) = 1.0 • 10 -3 . In our analysis we assign a 
30% error to this parameter so that the adopted range for L 5 reaches at the upper end the 
value implied by the p 4 fit and at the lower end the value obtained in the p 6 fit of Ref . [IE] ■ 
The combination (2L§ — L5) 00 can also be determined through resonance exchange. The 
only non-zero contribution comes from the exchange of pseudoscalar resonances. Within 
the SRA one gets 



(2L 8 -L 5 )°° = — ^pr « ^— = = —L? = -0.25 • 10" 3 . (5.13) 

8Ml 16MJ 4 ' 



The factor £?(yU SD ) in (|5.4j) and the 0(p 4 ) corrections A c , Aq' and A C A$ introduce 
additional dependences on the strong chiral couplings L 4 , L 6 and (3L 7 + L 8 ). At large -/V c , 



Lf = = and 



(4M| - 3M„ 2 - Ml)Fl 1 

v n V (it T ^ 00 nit: 1 n~ i 



(3^ + L *)°° = - u{M*-M*)* " 4 (2L§ " Ur = _0 ' 15 ' 10 ' (5 ' M) 

The same numerical estimate is obtained within the SRA, taking for the known con- 
tribution from rji exchange [THj . 



5.4 Strong couplings of 0(p 6 ) 

A systematic analysis of the LECs of 0(p 6 ) is still missing. Resonance contributions to 
some of the Xi have been studied in Refs. |1B1 Ej . 

Resonance dominance (that can be justified within large-iV c QCD) implies that the 
LECs of 0(p 6 ) occurring in the bosonization of the penguin operator Q e are determined 
by scalar exchange. The mass splitting in the lightest scalar nonet strongly influences those 
LECs. 
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We have estimated the relevant Xi with the scalar resonance Lagrangian discussed 
in Ref. 03] (setting gf = 0). The relevant resonance parameters in the nonet limit are 
Cd, c m , Ms, and e m , the latter governing the mass splitting within the scalar nonet. We use 
c m = Cd = F n /2, as determined from short- distance constraints |13j, and 

M s = 1.48 GeV, e s m = 0.2 (5.15) 

from a phenomenological analysis of mass spectra (these numbers correspond to sce- 
nario A of Ref. [IS!)- Even within resonance saturation, this is not a complete calculation 
of the relevant LECs of 0(p 6 ) but we expect it to capture the most significant physics. We 
refrain from reporting explicit numerics for the individual LECs here. Numerical values 
for the relevant combinations are reported in the next section. 

Finally, one can include nonet breaking effects within the framework of Ref. 03- m 
the chiral resonance Lagrangian, these effects are needed in order to understand the scalar 
mass spectrum (the coupling k m and 75 of Ref. [IH]). Once the resonances are integrated 
out, nonet breaking effects, sub-leading in 1/N C , appear in the Xi and therefore in the weak 
LECs g$Ni. Although this is far from being a complete analysis of sub-leading corrections it 
gives already an indication of their size. For all the quantities of physical interest, inclusion 
of k m and 75 produces shifts within our estimate of 1/N C corrections based on varying the 
chiral renormalization scale (see discussion above). 

5.5 Strong couplings of 0(e 2 p 2 ) 

Four combinations of the appear directly in the local amplitudes AqA^ of 0(e 2 G$p 2 ): 

Ui = K x + K 2 

U 2 = K 5 + K 6 

U 3 = K 4 -2K 3 

U 4 = K 9 + K 10 . (5.16) 

Within our large- iV c estimates, also other combinations of Ki appear through the couplings 
g89ewk and g&Zi. The ones relevant for K — > 1m decays are K 7 , K 8 , Kg, K w , Ku, Kyi, Ki3- 
It turns out that all the relevant combinations can be obtained from existing estimates 
|48| I49j. which we now briefly review. 

The LECs Ki can be expressed as convolutions of a QCD correlation function with the 
electromagnetic propagator. Therefore, their calculation involves an integration over the 
internal momenta of the virtual photon, which makes reliable numerical estimates difficult 
even at large iV c . In contrast to the strong LECs L\, the dependence of the Ul on the 
CHPT renormalization scale v x is already present at leading order in 1/N C . In addition, 
the Ki depend also on the short-distance QCD renormalization scale /xsd and on the gauge 
parameter £. Whenever numerical estimates are reported in the following, they refer to 
the Feynman gauge (£ = 1) and fisu = lGeV. 
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A first attempt to estimate the couplings Ki, using the extended Nambu-Jona-Lasinio 
model at long distances, has found the results 



[3 U[ + U r 2 ] 



K 



lOV^X 



M p 
M p 
M p 
M n 



(2.85 ±2.50) ■ KT 3 
-(2.5 ±1.0) ■ 10~ 3 
(2.7 ± 1.0) • 10~ 3 
(4.0 ± 1.5) ■ 10" 3 . 



The last two equations imply (adding the errors linearly): 



M p ) = -(1.3 ±2.5) ■ 10~ 3 



Moreover, in the limit N r 



oo, one has the relation [IB] 

Ul = 2U{ 

and the couplings K 7 ,K 8 are subleading. We therefore take K 78 (M p ) = 0. 
The remaining couplings needed were obtained at large N c in Ref. 



(5.17) 



(5.18) 



(5.19) 



through the 

evaluation of the relevant correlation functions in terms of narrow hadronic resonances. 
Within the SRA, one gets 19 : 

+ u 



K r 



-"■12 



-"■13 




, 27 33 , 

-i + — ln2 

^ 4 2 



In 2 



(5.20) 



Taking /x SD = 1 GeV, v x = M v and ( = 1, this gives K r xl = 1.26 • 10" 3 , K{ 2 = -4.2 • 10' 
and K[ 3 = 4.7 • 10 -3 . Inserting the SRA prediction from (pT27Ijl into (|P7j) . we get: 



U{{y x = M p ) = -5.0-10" 



Uliy x 



17.9-10" 



(5.21) 



A direct evaluation of U[ and L/J is in principle possible within the SRA ■ However, it 
requires an analysis of resonance couplings beyond the known results of Ref. [IE] . 



6 Numerical results 

We are now in the position to quantify the size of NLO contributions to the relevant isospin 
amplitudes, due to both chiral loops and local couplings in the effective theory. The master 
formulas for the amplitudes at NLO are given in Eq. ()4.6[) and (|4.7|) . They are organized 
in such a way as to easily identify the distinct sources of IB and to separate the LO from 
the NLO contributions in the chiral expansion. In Tables ^ 121 and [3] we report explicit 
results for the isospin amplitudes A n , n = 1/2,3/2,5/2, quoting for each component the 
following quantities: 
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• The LO contributions af^\ 

• The NLO loop corrections AlA)^\ consisting of absorptive and dispersive compo- 
nents. The dispersive component depends on the chiral renormalization scale v x 
(fixed at 0.77 GeV). 

• The NLO local contributions to the CP-even and CP-odd amplitudes, denoted re- 
spectively by [A C AW] + and [A C A^]~. Our estimates of [AcA^^ at the scale 
v x = 0.77 GeV are based on the leading 1/N C approximation. We discuss below the 
uncertainty associated with this method. 



The uncertainty in [AcA.^] 1 * 1 has two sources, related to the procedure used to estimate 
the NLO local couplings (see Sec. EJ), and we quote them separately in the tables. The 
first one corresponds to the short- distance input, essentially the renormalization scale used 
to evaluate the Wilson coefficients. We estimate this uncertainty by varying the scale fisD 
between 0.77 GeV and 1.3 GeV. The second one derives from working at leading order 
in the large-iV c expansion. At this order, there is a matching ambiguity because we do 
not know at which value of the chiral scale the estimates apply. Therefore, we vary the 
chiral renormalization scale v x between 0.6 and 1 GeV. The results show that the second 
uncertainty (long-distance) dominates over the first one (short-distance) in most cases. 
Moreover, one should keep in mind that the errors quoted for the [A^^A^] 111 are strongly 
correlated. In phenomeno logical applications we shall take such correlations into account. 

Some comments on the numerical results are now in order. From chiral power counting, 
the expected size of NLO corrections is at the level of ~ 0.2 — 0.3, reflecting M^/(ATrF n ) 2 ~ 
0.2. This estimate sets the reference scale in the following discussion. 



The definition of [AcA™]* is: 




(6.1) 




(6.2) 
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Table 1: Numerics for A1/2' a[ i 7 Al^J, Ac A 



( x ) 


"1/2 




\AcA\^\ + 


[AcA\^\ 


(27) 


9 


1.02 + 0.47 i 


0.01 ± ± 0.60 


0.01 ± ± 0.60 


(8) 


v/2 


0.27 + 0.47 i 


0.03 ± 0.01 ± 0.05 


0.17 ± 0.01 ± 0.05 


(e) 


2^ 
3^ 


0.26 + 0.47 i 


- 0.17 ± 0.03 ± 0.05 


1.56 ± 0.06 ± 0.05 


(7) 




- 1.38 


- 0.30 ± 0.05 ± 0.30 


- 12.6 ± 2.5 ± 0.30 


(Z) 


3 


- 1.06 + 0.79 i 


- 0.08 ± 0.01 ± 0.18 


0.17 ± 0.01 ± 0.18 


(g) 


2^ 
3 


0.27 + 0.47 i 


- 0.15 ± ± 0.05 


- 0.15 ± ± 0.05 



Table 2: Numerics for A 3/2 : aW, A L Aty, A c A : 



(X) 


a (X) 
"3/2 




[Ac4?] + 




(27) 


10 
9 


- 0.04 - 0.21 i 


0.01 ± ± 0.05 


0.01 ± ± 0.05 


(e) 


4 
3^ 


- 0.69 - 0.21 i 


- 0.15 ± 0.02 ± 0.50 


1.74 ± 0.06 ± 0.50 


(7) 




- 0.47 


0.59 ± 0.02 ± 0.10 


1.70 ± 0.35 ± 0.10 


(Z) 


4 
3 


- 0.86 - 0.78 i 


0.02 ± 0.01 ± 0.30 


0.16 ± 0.01 ± 0.30 


(g) 


2 
3 


- 0.50 - 0.21 i 


- 0.15 ± ± 0.20 


- 0.15 ± ± 0.20 



Table 3: Numerics for Az/2- a 5 / 2 , A L A 5 , 2 , Ac A, 



(X) 


a {X) 
"5/2 




\^cA^ 2 Y 


[Ac4?r 


(7) 




- 0.51 


- 0.20 ± ± 0.10 


- 0.11 ± 0.01 ± 0.10 


(Z) 




- 0.93 - 1.15 i 


- 0.14 ± 0.01 ± 0.40 


0.01 ± 0.01 ± 0.40 



The following pattern seems to emerge from our results. On one hand, whenever the 
absorptive loop correction is small, the dispersive correction is dominated by the local 
contribution. Therefore, it is rather sensitive to the chiral renormalization scale and to 
the values of LECs. In these cases, the size of NLO corrections is rather uncertain, at 
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least within the approach we follow here in evaluating the LECs. Extreme examples of 
this behaviour are provided by AcAfJ 2 (of little phenomenological impact) and AcAy 2 
(which is instead quite relevant phenomeno logically) . 

On the other hand, whenever the absorptive loop correction is large, the dispersive 
component is dominated by the non-polynomial part of the loops and it is relatively insen- 
sitive to the chiral renormalization scale and to the values of LECs. In all relevant cases 
we have checked that the absorptive component of AlA^ is consistent with perturbative 
unitarity. Therefore we conclude that in these cases the size of NLO corrections is rather 
well understood, being determined by the physics of final state interactions. Typical exam- 
ples of this behaviour are given by A L A^l 3 / 2 , which have an important phenomenological 
impact. 

We conclude this section with some remarks on apparently anomalous results. 

• A L A\ 

/2,3/2 is 0(1). As discussed above, the physics underlying this result is well 
understood, being related to the absorptive cut in the amplitude. The key point is 
that this feature is absent at LO in the chiral expansion. It first shows up at NLO, 
setting the natural size of the loop corrections. NNLO terms in the chiral expansion 
are then expected to scale as NLO x(0.2 — 0.3), since corrections to the absorptive 
cut behave this way. Therefore A L Ayl 3 / 2 ~ 0(1) does not imply a breakdown of 
the chiral expansion. 

• [a c a ( ; 

72,3/2] * s 0(1). This result is determined by the large numerical coefficients 
multiplying the couplings N£ 78l3 , which turn out to have natural size within the 
leading 1/N C approximation. We observe, however, that in the case of the AI = 3/2 
amplitude (phenomenologically relevant) the leading 1/N C approximation is afflicted 
by a large uncertainty due to high sensitivity to u x . This uncertainty mitigates the 
apparent breakdown of chiral power counting. 

• Yet another surprising result is the one for [Ac^A^] - . The underlying reason is in 
the large size of the CP violating component of the Wilson coefficients O 9i i . Again, 
the operators 0,9,10 only make their first contribution at NLO in the chiral expansion. 



7 Phenomenology I: CP conserving amplitudes 

This section is devoted to a phenomenological analysis of K — > im decays including all 
sources of isospin breaking. The theoretical parametrization of the amplitudes is based 
on the NLO CHPT analysis discussed in the previous sections. Our goal is to extract 
information on the pure weak amplitudes (or equivalently on the couplings g$ and (727) and 
to clarify the role of isospin breaking in the observed K — > irir rescattering phases. All 
along we keep track of both experimental errors and the theoretical uncertainties related 
to our estimates of the NLO couplings at leading order in the 1/N C expansion. 
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7.1 Including the radiative modes 



When considering electromagnetic effects at first order in a, only an inclusive sum of 
K — ► 7T7T and K — ► 7r7r7 widths is theoretically meaningful (free of IR divergences) and 
experimentally observable. 

We denote the appropriate observable widths by T n ^(uj) for n = H — , 00, +0. These 
widths depend in general on the amount of radiative events included in the data sample, 
according to specific experimental cuts on the radiative mode. This dependence is com- 
pactly represented by the parameter uj. Denoting by A n the IR finite amplitudes as defined 
in Eq. ()2.3j) . the relevant decay rates can be written as 

r n[7] H = \A n \ 2 $„ G n (u) . (7.1) 



Here ^fs^ is the total cms energy (the appropriate kaon mass) and $ ra is the appropriate 
two-body phase space. The infrared factors G n (u>) are defined as 



G n (u) = 1 + - 

7T 



2tt Re J B„(M 7 ) + I n (M 7 ; u) 



(7.2) 



Note that G n (uj) is different from 1 only in the K° — > 7r + 7r~ and fC + — > 7T + 7r° modes. The 
factor £? n (M 7 ) arises from the IR divergent loop amplitude (its definition for n = H — is 
given in Eq. (|4.5j) ). while a/ir J n (M 7 ; a;) is the X — > tttt^ decay rate normalized to the 
non-radiative rate. The latter term depends on the treatment of real photons (hence on 
uj) and is infrared divergent. The combination of IR divergences induced by virtual and 
real photons cancels in the sum, leaving the uj dependent factor G n (uj). 

We discuss here in some detail the expression for G + -(uj), which plays an important 
phenomenological role. On the other hand, the inclusion of G + o{uj) only produces an effect 
of order (or aG^) and therefore represents a sub-leading correction. Its numerical 

effect will be taken into account, following the analysis of Ref. The explicit form of 
B_i__(M 7 ), the virtual photon contribution to G+-(uj), can be found in Eq. (|B.1|) . The real 
photon contribution a/ft J n (M 7 ; uj) arises from the decay K°(P) — > 7r + (p + )ir~ (p-)^(k) and 
it has the form 

J + _(M 7 ;^) = / cfe/ + _( S ;M 7 ) (7.3) 

M K\ 1 MT 



where 



s = (p++p_r, s min = 4M 2 , s max = (Mr- — M 7 ) (7.4) 

X ± = \(M* K -8- M 2 ) ± - ^AV 2 (M|, Sj M 2 ) (7.6) 
A(x, z) = x 2 + y 2 + z 2 — 2{xy + + 1/2) . (7.7) 
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The infrared divergence comes from the upper end of the integration in the dipion invariant 
mass (s ~ s max ). We have verified by analytic integration in the range M K [M K — 2u) < 
s < s max that J + _(M 7 ; u) has the correct M 7 dependence to cancel the infrared singularity 
generated by virtual photons. For uj/Mk <C 1, the analytic expression of 7 + _(M 7 ;u/) can 
be found in Eq. (21) of Ref. 0. The corresponding function G + -{uj) is plotted in Fig. 2 
of Ref. 0. 

Recently, the KLOE collaboration has reported a high-precision measurement of the 
ratio r + _/T o P, where the result refers to the fully inclusive treatment of radiative 
events. In order to use the KLOE measurement in our analysis, we need to calculate 
the fully inclusive rate (no cuts on the tnr'y final state). We have done this by numerical 
integration of Eq. (|7.3J) and we find 

1 + 0.67- 1(T 2 . (7.8) 



inclusive 



7.2 Constraints from measured branching ratios 

CP conserving K —>■ tttt phenomenology is based on the following input from Eq. (j7.1|) for 
n = +-,00,+0: 

'-M^f (7 - 9) 

It is convenient to express these equations in terms of the isospin amplitudes A ,A 2 , A 2 
and the phase shift Xo ~ X2 (as defined in Eq. (j2.2|) V With r = (C + _/C o) 2 one obtains 5 

A\ = g C+o 

(A ) 2 + (A 2 ) 2 = 2 -Cl_+ l -Cl (7.10) 

A2 , , r-l + (^) 2 (2r - ±) 

— cos(xo - X2) 



A VAU AZ/ V2(l + 2r) 

In general, in the presence of isospin breaking, these three independent experimental con- 
straints are not sufficient to fix the three isospin amplitudes (A , A 2 , A 2 or A\/ 2 , ^.3/2, ^4.5/2) 
plus the phase difference (xo — X2)- In the previous sections, we have shown how CHPT 
relates the amplitude ^4 5 / 2 to A\/ 2) thus effectively reducing the number of independent 
amplitudes. Including also all other isospin breaking effects, we can extract the couplings 
g 8 and g 27 from Eqs. (|7.1()jl. 



7.3 CHPT fit to K -> 7T7T data 

Using Eqs. (I7.10J) as starting point, we perform a fit to g 8 , g 27 and the phase difference 
Xo — X2- In order to do so, we employ the CHPT parametrization for the Aj. The detailed 

5 Note that in the last equation one has to use Ai and not (as done in the isospin conserving 
analyses). For this reason the extraction of the phase shift is related to the A J = 5/2 amplitude. 
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relations between ^.1/2,3/2,5/2 (presented in Sees. El EJ and ^0,^2,^2", to first order in 
isospin breaking, are reported in Section 18.11 We leave the phase difference Xo ~ X2 as a 
free parameter because one-loop CHPT fails in reproducing the strong s-wave phase shifts. 

Apart from g 8 and g 27 , the amplitudes A n depend on the LO couplings (fe^cwk, Z and 
on a large set of NLO couplings. Given our large-iV c estimates for (fefiWk and for the 
ratios of NLO over LO couplings (gsNi)/g 8 , ■■■ (see Sec. EJ), we study the constraints 
imposed on gg,g 27 by the experimental branching ratios. In this process we keep track of 
the theoretical uncertainty induced by the use of a specific approximation in estimating 
the relevant LECs (leading order in the large- expansion). In practice, this reduces to 
studying the dependence of the amplitude (and of the output values of g 8 and #27) on two 
parameters: the short distance scale (/^sd) and the chiral renormalization scale (v x )- 

In summary, the experimental input to the fit is given by the three partial widths 
r + _ ,00, +0 (kaon lifetimes and branching ratios) [TO] and by the new KLOE measurement 
for r + _/T o [§]■ The theoretical input is given by the NLO CHPT amplitudes as well as the 
estimates for g ewk and the NLO couplings. As primary output we report Re g 8 , Re g 27 and 
Xo ~ X2- Derived quantities of interest for phenomenological applications will be reported 
subsequently. 

(1) Using the NLO isospin conserving amplitudes (IC-fit), we find 

Reg 8 = 3.665 ± 0.007 (exp) ± 0.001 (/i SD ) ± 0.137 (u x ) 
Reg 27 = 0.297 ±0.001 (exp) ± 0.014 (u x ) (7.11) 

X0-X2 = 48.6 ±2.6 (exp) . 

Using instead the tree-level (LO) amplitudes in the isospin limit would lead to Reg% = 
5.09±0.01 and Re #27 — 0.294±0.001. This result is in qualitative agreement with Ref. [50] : 
NLO chiral corrections enhance the 1 = 1/2 amplitude by roughly 30%. 

(2) Using the full NLO isospin breaking amplitudes (IB-fit), we find 6 

Re £ 8 = 3.650 ±0.007 (exp) ±0.001 (/i SD ) ± 0.143 (v x ) 
Reg 27 = 0.303 ±0.001 (exp) ± 0.001 (/i SD ) ± 0.014 (v x ) (7.12) 
X0-X2 = 54.6 ±2.2 (exp) ±0.9 (u x ) . 

In the IB case, a tree-level (LO) fit leads to Re# 8 = 5.11 ±0.01 and Reg 27 = 0.270 ±0.001. 
A few remarks are in order: 

• Using NLO amplitudes, both g$ and g 27 receive small shifts after inclusion of IB 
corrections. While this could be expected for g 8 , it results from a cancellation of 
different effects in the case of g 27 (at tree level the inclusion of isospin breaking 

6 The uncertainty in = (1.061 ± 0.083) • 10~ 2 produces errors one order of magnitude smaller 
than the smallest uncertainty quoted above. 
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reduces g 2 7 by roughly 10% ). Note also that competing loop effects reduce the v x - 
dependence of Re g 2 j (IB-fit) to only ±0.002. As a more realistic estimate of the 
long-distance error we have chosen to quote the ^-dependence induced by each one 
of the competing effects (for example the isospin-conserving loops). 

The output for g 8 and g 27 is sensitive to the input used for the strong LECs Lj of 
0{j) A ). The results of Eqs. (|7.11|) and (j7.12j) correspond to the central values quoted in 
Sec. |S] We have repeated the fit with non-central input and have found the variations 
in g 8 and (727 to be below 5%. 

In obtaining the results in Eqs. (|7.11j) and ()7.12|) we have used the large- N c predictions 
for the ratios {gsN^/gg. We have also employed the alternative procedure of using 
large N c directly for the couplings gsNi. In this case we find g s = 3.99 and g 27 = 
0.289, reflecting the change in size of the p A local amplitudes. All other quantities of 
phenomenological interest are stable under this change in the fitting procedure. 

Some derived quantities of phenomenological interest are the ratios of isospin ampli- 
tudes: ReA /ReA 2 , ReAg/ReA^ , f 5 / 2 = ReA 2 /ReA% — 1. From our fit we find 



Rev4 

ReA 2 
ReA 



ReA 2 J IB-fit 



IB -fit 



20.33 ±0.07 (exp) ±0.01(// S d) ±0.47(z/ x ) 

22.09 ±0.09 (exp) ±0.01(// S d) ±0.05(z/ x ) (7.13) 



/ 5 / 2ffi _ fit = (8.6 ±0.03(exp) ±0.01(/i SD ) ± 2.5 (zaJ) • 10 



-2 



In the absence of isospin breaking, one finds instead f^/ 2 = and ReAo/Re/^ 
22.16 ±0.09. 



7.4 Isospin breaking in the phases 

This section is devoted to understanding isospin breaking in the rescattering phases of 
K — > 7T7r. If isospin is conserved Watson's theorem predicts Xo ~ X2 = — ^2 ~ 45°. For 
a long time this prediction has not been fulfilled by the data, as one typically encountered 
Xo - X2 ~ 60°. 

The situation has recently improved. Using the KLOE data jU] and working in the 
isospin limit, our fit (|7.11|) gives xo ~ X2 ~ 49° and so there seems to be no more phase 
problem. However, the inclusion of isospin breaking appears to reintroduce the issue. In 
order to understand what is going on, we analyse in detail the various factors determining 
Xo - X2- 

The last of Eqs. ()7.10|) can be rewritten as 

At r-l±(^) 2 (2r-i) 

^ cos(xo - x2H v2(i + 2;)(i +/5/ 2) • (7 - i4) 
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Let us first calculate the right-hand side without an I = 5/2 amplitude (f$/2 = 0). With 
the old value r = 1.1085 as input (PDG2000 [52 ), one obtains 

4 L cos(xo-X 2 )= 0.02461. (7.15) 

With A? /A = 0.045 (based on the IC-fit), this leads to the standard puzzle that Xo ~ X2 
is much bigger than 45°: 

Xo-X2 = 57°. (7.16) 

What could be the reasons for this discrepancy of about 30 % (cos 45°/ cos 57° = 1.30)? 
Let us consider several effects: 

• First of all, the right-hand side of ()7.14j) has changed with the recent KLOE result 
[0] r = 1.1345 to give 

4^ cos( X o-X 2 )= 0.02987. (7.17) 

This is a sizable correction of about 18 % and it goes more than half-way in the right 
direction to decrease the phase difference. 

• Taking into account isospin breaking introduces an / = 5/2 amplitude via the ratio 
f 5 / 2 in Eq. ()7.14|) . According to our results ()7.13|) . 

f 5/2 = (6.5 (loops) + 2.1 (local) ± 2.5(z/ x )) • 10~ 2 , (7.18) 

increasing again the discrepancy. This value is dominated by loop contributions. 
Even if one changes the sign of the relevant combination of LECs, f 5 / 2 would still 
be positive. Note that ()7.18|) amounts to a correction of ~ 8 % (in the "wrong" 
direction). As already noted by Wolfe and Maltman jjj, it seems impossible to solve 
the phase problem with a reasonable choice of counterterms. 

• Finally, the "infrared factor" for the H — mode must be taken into account. This is 
straightforward with the inclusive measurement of KLOE. We find r = 1.127, which 
increases again the discrepancy. Including also the effect of f$/2, we obtain 

4^ cos( X o-X 2 )= 0.02611, (7.19) 
leading to X o - X2 = (54.6 ± 2.4)°. 

Before addressing the question whether this result is in disagreement with the nir phase 
shift prediction [53J 5$ — 62 = (47.7 ± 1.5)°, it is mandatory to study the effect of isospin 
breaking on the phases themselves. We use the general decomposition 

Xi = Sj + 7/ (1 = 0,2) , (7.20) 

where 7/ represents an isospin breaking correction. The 7/ are related to isospin breaking 
dynamics in irn rescattering as well as to the presence of radiative channels IM] • Since 
the analysis of Ref. [Sj, new information on radiative corrections in tttt scattering has 
become available, allowing for a reevaluation of 70 — 72- 
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7.5 Optical theorem and 70 — 72 



The K — > 7T7T amplitudes at NLO in CHPT allow for a perturbative evaluation of 70,2- We 
find 7 

70 = (-0.18 ±0.02)° 

7 2 = (3.0 ±0.4)°, (7.21) 

where the error is obtained by varying the chiral renormalization scale v x as in the main 
fit. Setting the NLO local terms to zero would lead to results within the range quoted in 
Eq. (|7.2ip . This evaluation incorporates the constraints of the optical theorem at leading 
order in perturbation theory. In practice, this only reflects the O(e 2 p ) mixing between the 
1 = and I = 2 tttc channels, completely missing both higher-order corrections and the 
new physical effect due to the radiative channel irnj. In order to improve upon these per- 
turbative results, a more general analysis of the optical theorem for K° — > tttt amplitudes is 
required. We shall follow here the approach of Ref. J5J, except for a few details. The main 
novelty lies in the final stage, in which one needs an explicit calculation of isospin breaking 
effects in im scattering: we use the results obtained at 0(e 2 p 2 ) in CHPT in Refs. 155]. 

We now summarize the steps involved in the optical theorem analysis of Ref. [3] , rele- 
gating some technical details to App. |Dj For this section, CP is assumed to be conserved. 

1. The first step is to work out the consequences of the optical theorem for K° — ► ttti 
amplitudes, considering the following intermediate states: 7r + 7r~, 7r°7r° and 7r + 7r~7. 
For the radiative amplitudes describing K° — ► 7r + 7r~7 and 7r + 7r _ 7 — > tttt we use the 
leading Low parametrization, thus neglecting possible structure dependent terms 8 . In 
this approximation the radiative amplitudes are known in terms of the non-radiative 
ones. Under the assumptions listed above, and collecting the K° — > tttt amplitudes 
in a two-component vector A, the optical theorem has the following form: 

AbsA = ft (T ] + U) A (7.22) 

where ft = \Jl — 4M^ /M|-, while T and 1Z are two-by-two matrices: T is related to 
the s-wave projection of the tttt T-matrix, while TZ encodes the effect of both radiative 
modes and phase space corrections induced by mass splitting. 

The explicit form of Eq. ()7.22|) is best derived by working with tttt states in the charge 
basis (7r + 7r~, 7r°7T ) where it is more transparent to deal properly with IR singularities 
and phase space corrections. Special care is needed in removing the IR and Coulomb 
singularities from the amplitudes A+-, 7+-, 00 an d Tf_ )+ _. This step involves an 
arbitrary choice, which only affects the intermediate states of the analysis but not 

7 The results depend on the ratio 58/327, for which we use the IB-fit output. 
8 This is known to be an excellent approximation for K — > n + ir~^f. 
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the final results. We adopt the following prescription 9 : 



A(K° - 


• 7T + 7T~) 


= A+- 


A(n°n° - 


' 7T + 7T~) 


= ^+-,00 


A(n + 7i~ — 


> 7T + 7r~) 





exp {a -8^} 
exp{a5^r} 



(7.23) 



where the infrared singularity is separated in the factors B n7T and C nn , whose form 
is reported in App. [01 These factors depend only on the charges and kinematical 
configuration of the external particles. 

2. In order to make contact with standard treatments, it is convenient to represent 
Eq. (J7.22JI in the "isospin" basis for tttt amplitudes. Explicit relations between charge 
and isospin amplitudes are reported in Ref . jS] . In the isospin basis the matrices have 
the form 



T 



1 T T ^ 

-J 00 ^02 



T~2Q ^22 



K 



( ^(A + _T *+^ 



^ Y (A+-T 2 * + 5 A 



^(A+_T 00 + 5 H 



(a+_t 2 * + 5, 



J 

(7.24) 



where, using the notation 



(f) 



d(cos6) f(cosO) 



the various quantities have the following structure: 

1 



ab 



A 4 

6 4 



64tt 

mi 



+ 2a Re (B ww ) + 



K 



32tt 



fx 



rad 
2 



rad 



(7.25) 



(7.26) 



Xb is the s-wave projection of the b — > a tttt scattering amplitude. The factor A + _ 
receives a contribution from phase space corrections (pion mass splitting), one from 
virtual photons and one from real photons (/™ d ). Likewise, £+_ reflects both 

virtual corrections (C^) and real photon effects (/ 2 ad ). The definition of the phase 
space factors and d$ + _ 7 , as well as of /{ ad is reported in App. |D| We remark 

here that A + _ and <5+_ are free of infrared singularities, as these cancel in the sum 
of real and virtual photon contributions. 



9 The 7T7T amplitudes are functions of two of the three Mandelstam variables s,t,u. In the following 
we set s — Mjj and trade the other independent variable for the cms scattering angle 6. Moreover, the 
explicit dependence on cos# is suppressed in order to keep the expressions compact. 
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3. At this point one needs a general parametrization of the matrix 7}j, the T- matrix 
restricted to the dimension-two subspace of tttt channels. Assuming T-invariance 
(but not unitarity of the S-matrix restricted to this subspace), an explicit form is 
given by 



T 



1 



\ 



2i 



a e 



{rj 2 e 



2i& 2 



r 



(7.27) 



2i 



J 



in terms of five parameters (two phase shifts, two inelasticities and one off-diagonal 
amplitude). If one assumes that only one extra state couples to the ones considered 
here (namely the 7r + 7r~7 state), then the inelasticities are correlated, as noted in 
Ref. [54]. Since our subsequent discussion does not depend on the inelasticities, we 
do not elaborate further on this point. 

4. The next step is to assume an ansatz for the K — > tttt amplitudes of the type 

Aj = A/ e 4{5 ' +7l) (/ = 0, 2) (7.28) 

and work out the constraints imposed upon 77 by Eq. (|7.22j) . Solving for sin 70 and 
sin 72 to first order in a and taking into account A 2 /A -C 1, one finds [5] 



sm7o 
sin 72 



(3 (Re(7^oo) ~ tan<5 Im(7^ o)) — O(asm8 ) 
1 



<5 



T20I + 



COS(^ 



(Re(72-2o) cos^o ~ Im(72-2o) sin5 



(7.29) 



The key feature of Eq. ()7.29j) is that in the expression for 72 the isospin breaking 
effects get once again enhanced by the factor A /A 2 ~ 22. 

5. The final step consists in evaluating T 02 , A + _ and for which we need an explicit 
expression for the tttt amplitudes with isospin breaking |B31 15TI] . as well as explicit 
expressions of B wn and C n7T . The details of the calculations cannot be given in a 
concise way and we report here only the results. 
For 7q2 we find 

V2(M% ± -M* ) 



Using the results of Refs. jSHl EE] and their estimate of the relevant LECs Ki, we 
obtain 

Ao2 p2 = (0.78 ± 0.83) + i 0.54 (7.31) 
where we have added the various errors in quadrature. 

For the radiative factors the calculation cannot be done in a fully analytic form and 
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we employ Monte Carlo integration to deal with the real-photon contribution to <5 + _ . 
We find 

A + _ = -0.81 ■ KT 2 

S + _ = 0.09 • 10~ 2 . (7.32) 
Using these input values in Eq. ([7.29)1 . we arrive at 

7o = -0.2° 

7 2 = (6 ± 3)° . (7.33) 

The conclusion is that the optical theorem estimate of 70 — 72 agrees roughly with 
the perturbative estimate (j7.21|) and that it tends to worsen the discrepancy between the 
theoretical prediction of 5q — 62 [SB] and its phenomenological extraction from K — ► ttti 
decays. Explicitly one has 

(S Q -S 2 )^ n = (47.7 ± 1.5)° 

(S - 8 2 ) K ^ = (60.8 ± 2.2 (exp) ± 0.9 (v x ) ± 3.0 (72))° • (7.34) 

Although the precise KLOE measurement of the ratio of K$ — > 7T7t rates has considerably 
improved the situation we still obtain a difference of about 13° for the phase shift difference 
in the isospin limit between the two determinations in Eq. ([7.34)1 . The theoretical error 
is much bigger in the present case due to uncertainties in the NLO LECs. However, we 
observe that more than half of this difference is due to the AI = 5/2 loop amplitude that 
depends only on the well-established lowest-order electromagnetic LEC Z in the Lagrangian 
([2.8)1 . In order to obtain a phase shift difference in the isospin limit below 50°, the local 
amplitude with AI = 5/2 would have to be more than twice as big and of opposite sign. 
While such an explanation cannot be totally excluded at this time, the discrepancy in the 
two entries of Eq. ([7.34)1 certainly warrants further study. 

The AI = 5/2 amplitude induced by isospin violation in the octet amplitude is small 
because it only arises at NLO and it is of purely electromagnetic origin. One may wonder 
whether isospin violation in the 27-plet amplitude, which occurs already at leading order, 
could compete. Whereas isospin violating contributions to the AI = 1/2,3/2 amplitudes 
proportional to G27 are certainly negligible, the effect on the AI = 5/2 amplitude is worth 
investigating. 

It is straightforward to calculate isospin violation in the LO 27-plet amplitudes in ([3.6)1 
due to mass differences and 7r° — r\ mixing. We are only interested in the resulting AI = 5/2 
amplitude, entirely due to the quark mass difference: 

^—G^Ml-MDe^ . (7.35) 
This amplitude may now be compared to the corresponding 5/2 amplitude in ([4.6)1 : 

A/2 = -e 2 G 8 F^ (M + ZAfh . (7.36) 
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With the numerical information of Table El and Eq. ([7.12)1 . we obtain for the ratio 

A {27) 

■" 2 ~ 6 • 1(T 2 . (7.37) 



Visp A 5 / 2 

The conclusion is that the 27-plet contribution to the AI = 5/2 amplitude is of the same 
sign and only about 6 % of the octet contribution. The impact on the AI = 1/2,3/2 
amplitudes is of course much smaller still. Isospin violation in the 27-plet amplitude can 
safely be neglected. 



8 Phenomenology II: CP violation 

The main contents of this section have already been published in Ref. jHZ|. They are 
included here for completeness. 

8.1 Isospin violation and e' 

The direct CP violation parameter e' is given by 

Tim4 Iim4 2 



e ' = L e i(x2-xo) Rey42 



ReA ReA 2 



•1) 



y/2 Rev4 

The expression (18.1)1 is valid to first order in CP violation. Since ImAi is CP odd the 
quantities ReAj and xi are om y needed in the CP limit (/ = 0, 2). 

To isolate the isospin breaking corrections in e', we write the amplitudes Aq, A2 more 
explicitly as 

A 2 e^ = 4/2 + ^3/2+^5/2, (8.2) 

where the superscript (0) denotes the isospin limit and SAi/2,3/2, ^5/2 are first order in 
isospin violation. In the limit of isospin conservation, the amplitudes Aai would be gen- 
erated by the AI component of the electroweak effective Hamiltonian. 

To the order we are considering, the amplitudes Aai have both absorptive and disper- 
sive parts. To disentangle the (CP conserving) phases generated by the loop amplitudes 
from the CP violating phases of the various LECs, we express our results explicitly in terms 
of Visp Aai and Abs Aai- Writing Eq. (|8.2j) in the generic form 

Aje 1 ^ =A n = A^+6A n , (8.3) 

we obtain to first order in CP violation: 



ReAj = sj (Re[Visp A n }) 2 + (Re[Abs A n }) 2 

ImAj = (ReA^ 1 (\m[VispA n }Re[VispA n }+lm[AbsA n }Re[AbsA r ]) (8.4) 
e iXl = (ReA/) -1 (Re[DispAn]+iRe[AbsAn\) . 
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Using the second equality in Eq. (|8.3|) . one can now expand ReAi and ImA/ to first order 



in isospin breaking. With the notation \A^\ = \j (Re[Disp An^}) 2 + (Re[Abs An } }) 2 , we 



i (o)i 



find 

ReA 7 
Iim4 7 



^n } \ + \^n ] \~ (MVispA^] Re[Disp5An] + Re[AbsA^} Re[Abs5A n }) (8.5) 
^i 0) | _1 {impispA^} Re[DtspA^] + lm[AbsA^] Re[AbsA^}} 
+ |^i 0) | _1 {lm[Disp5An\ Re[VispA^\ + hn[VispAf] Re[Visp5A n ] 
+ lm[Abs 5 A n ] Re[Abs A®] + lm[Abs A$] Re[Abs 5A n ] } 

- |^i 0) f 3 [Repisp A®] Re[Visp 5A n ] + Re[Abs A^} Re[Abs 5A n ] } x 

[\m\VispAf] Re[VispA^} + lm[Abs A®] Re[Abs A ( ° ] }} , (8.6) 

where the first term in each equation above represents ReAf^ and lmAf\ respectively. 

We now turn to the different sources of isospin violation in the expression (j8.1|) for e'. 
We disregard the phase which can be obtained from the K — > branching ratios. The 
same branching ratios are usually employed to extract the ratio ujs = ReA 2 /ReAQ assuming 
isospin conservation. Accounting for isospin violation via the general parametrization (|2.3|) . 
one is then really evaluating uo + = ReA 2 /ReA rather than u>s- The two differ by a pure 
AI = 5/2 effect: 



fb/2 



UJ + (l + 

ReA 2 
ReAt ~ 



1.7) 



1 . 



Because u> + is directly related to branching ratios it proves useful to keep u + in the nor- 
malization of e', introducing the AI = 5/2 correction / 5 / 2 [5J. 

Since Iim4 2 is already first order in isospin breaking the formula for e' takes the following 
form, with all first-order isospin violating corrections made explicit: 



Ax2-xo) 



UJ+ 



ImA 



(o) 



ReA 



(0) 



5/2, 



lmA 2 
Re^F 



where 



(o) 



A 



lmA ReA y Q 
ImA^ ReA 



- 1 



19) 



(8.10) 



With the help of Eqs. ()8.5|) and (|8.6|) . one obtains 



1mA, 



ImA, 



(0) 



Af/ 2 \ 1 {lm[VispA^ 2 ]Re[VispA^ 2 }+lm[AbsA^ 2 }Re[AbsA^ 2 ]} (8.11) 



A 



(o) 

3/2 



-1 



{lm[Disp (<L4 3 /2 + A/2)] Re[VispA 



(0)1 
3/2J 
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+ lm[Abs (dA 3/2 + A 5/2 )] Re[AbsAf/ 2 }} 



[8.12) 



A 



(0) 
1/2 



+ 



Re[VispA^/ 2 ] Re[VispSA 1/2 }+Re[AbsAY/ 2 } Re[AbsSA 1/2 
[lm[Visp Afj 2 ] Re[Visp Af/ 2 ] + lm[Abs Af/ 2 ] Re[Abs Af/ 2 \] ~* x 
{lm[Visp5A 1/2 } Re\VispAfJ 2 ] + lm[Visp AfJ 2 ] Re[Disp 5A x/2 \ 
+ lm[Abs 5A 1/2 ] Re[Abs Af/ 2 ] + lm[Abs Af} 2 ] Re[Abs 5A 1/2 ] } 



/5/2 



A (0) 



-2 



(8.13) 

{Re[VispA { °/ 2 } Re[VispA 5/2 ] + Re[Abs A$} 2 ] Re[Abs A 5/2 }} . (8.14) 



These expressions are general results to first order in CP and isospin violation but they 
are independent of the chiral expansion. Working strictly to a specific chiral order, these 
formulas simplify considerably. We prefer to keep them in their general form but we will 
discuss later the numerical differences between the complete and the systematic chiral 
expressions. The differences are one indication for the importance of higher-order chiral 
corrections. 

Although Im^2 is itself first order in isospin breaking we now make the usual (but 
scheme dependent) separation of the electroweak penguin contribution to IH1A2 from the 
isospin breaking effects generated by other four-quark operators: 

ImA 2 = ImA e 2 mp + ImAr _cmp • (8.15) 

In order to perform the above separation within the CHPT approach, we need to identify 
the electroweak penguin contribution to a given low-energy coupling. In other words, we 
need a matching procedure between CHPT and the underlying theory of electroweak and 
strong interactions. Such a matching procedure is given here by working at leading order 



in 1/N C . Then, the electroweak LECs of 0(G 8 e 2 p r ' 



(n 



0,2) in ImAg 011 emp must be 



calculated by setting to zero the Wilson coefficients CV, C$, Cg, C10 of electroweak penguin 
operators. 

Splitting off the electromagnetic penguin contribution to ImA 2 in this way, we can now 
write e' in a more familiar way as 



where 



V2 



e i(X2- X o) u+ 



ImA ( 



(0) 



ReA( 0) 



1-fieff)- 



lmAl mp 



eft 



^IB — A — /5/2 



ReA 



(0) 



IB 



ReAf 



ImA° on - emp 
ImAp 



(8.16) 



(8.17) 



The quantity f2 e ff includes all effects to leading order in isospin breaking and it generalizes 
the more traditional parameter Although Q\b is in principle enhanced by the large 
ratio ReAp /ReA^ the actual numerical analysis shows all three terms in (|8.17j) to be 
relevant when both strong and electromagnetic isospin violation are included. 
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8.2 Numerical results 

We present numerical results for the following two cases: 

i. We calculate fl e s and its components for a = 0, i.e., we keep only terms proportional 
to the quark mass difference (strong isospin violation). In this case, there is a clean 
separation of isospin violating effects in Iul4 2 . We compare the lowest-order result 
of 0{m u — m d ) with the full result of 0[{m u — m d )p 2 }. 

ii. Here we include electromagnetic corrections, comparing again 0(m u — md, e 2 p°) with 
0[{m u — m ( i)p 2 ,e 2 p 2 }. In this case, the splitting between ImAg 11113 and ImA2° n ~ cmp is 
performed at leading order in 1/N C . 

The LO entries depend on Reg 8 /Reg27 as well as on Im(g 8 g e wk) /Im#8- Subleading effects 
in 1/N C are known to be sizable for the LECs of leading chiral order. We will therefore not 
use the large- iV c values for Reg 8 , Re#27 in the numerical analysis but instead determine 
these couplings from our fit to the K — > ttti branching ratios. 

The other combination of interest is the ratio Im(g 8 g ew ]J/lHU7 8 . In this case, existing 
calculations beyond factorization jSH] suggest that the size of 1/N C effects is moderate, 
roughly — (30 ± 15)%. As a consequence, it turns out that the dominant uncertainty comes 
from the input parameters in the factorized expressions. Finally, one also needs the ratio 
lm(g$g ew k) non ~ emp /Img$: in this case, leading large-iV c implies —3.1 ± 1.8 (error due to 
input parameters), while the calculation of Ref. PP gives —1.0 ± 0.5. Given the overlap 
between the two ranges and the large error in the large- N c result, we use in the numerics 
the range implied by leading large- N c . 

At NLO the quantities we need to evaluate depend on the ratio of next-to-leading to 
leading-order LECs. In Table EJ we use the leading l/N c estimates for the ratios GgNi/Gg, 
. . . The final error for each of the quantities fire, A , /s/2 and Q e ff is obtained by adding 
in quadrature the LO error and the one associated to weak LECs at NLO. Moreover, 
only / 5 / 2 and ReA^ /ReA^ depend on the ratio ggjgn- In these cases we rely on the 
phenomenological value implied by our fit. Some of the errors in Tableware manifestly 
correlated, e.g., in the LO column for a^O. 

The NLO results are obtained with the full expressions (|8.11jl . . . . , (|8.14jl . Using instead 
the simplified expressions corresponding to a fixed chiral order, the modified results are 
found to be well within the quoted error bars. We therefore expect our errors to account 
also for higher-order effects in the chiral expansion. 

We have also employed an alternative procedure for estimating the non-leading weak 
LECs. In contrast to the previous analysis, we now apply large N c directly to the LECs 
GsNi, . . . This amounts to assuming that the failure of large N c for Gg is specific to the 
leading chiral order and that the non-leading LECs are not significantly enhanced compared 
to the large-iV c predictions. Of course, this implies that the local amplitudes of OiGgp 4 ) 
are less important than in the previous case. Consequently, the fitted value for g s comes 
out quite a bit bigger than in Eq. (|7.12p . whereas #27 gets smaller (see Sec. 17. 3J) . However, 
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LO 


a = 

LO+NLO 


a 

LO 


7^0 

LO±NLO 




11.7 


15.9 ±4.5 


18.0 ±6.5 


22.7 ± 7.6 


A 


- 0.004 


- 0.41 ±0.05 


8.7 ±3.0 


8.4 ±3.6 


h/2 











8.3 ±2.4 


^eff 


11.7 


16.3 ±4.5 


9.3 ±5.8 


6.0 ±7.7 



Table 4: Isospin violating corrections for e' in units of 10~ 2 . The first two columns refer 
to strong isospin violation only (m u ^ rrid), the last two contain the complete results 
including electromagnetic corrections. LO and NLO denote leading and next-to-leading 
orders in CHPT. The small difference between the value of /s/2 reported here and the 
one in Eq. (|7.13|) is due to higher-order effects in isospin breaking (absent in this table 
according to Eq. (|8.14|0 . 

the isospin violating ratios in Table 0] are very insensitive to those changes. Not only are 
the numerical values in this case well within the errors displayed in Tabled] but they are 
in fact very close to the central values given there. 

Finally, we have investigated the impact of some sub leading effects in 1/N C [45J. Al- 
though this is not meant to be a systematic expansion in 1/N C , the nonet breaking terms 
considered in [321 may furnish another indication for the intrinsic uncertainties of some of 
the LECs. The size of those terms depends on the assignment of isosinglet scalar reso- 
nances. Since nonet breaking effects are large in the scalar sector they affect most of the 
entries in Table 0] in a non-negligible way, although always within the quoted uncertain- 
ties. Employing scenario A for the lightest scalar nonet 03], Q e s in ()8.16|) decreases from 
6.0 • 10~ 2 to -1.4 • 10~ 2 . 

The lessons to be drawn from our analysis of isospin violating corrections for e' are 
straightforward. Separate parts of those corrections turn out to be sizable. A well-known 
example are the contributions of strong isospin violation to 7r° — rj mixing where the sum 
of r) and rj' exchange generates an f2iB of the order of 25 % [HH] . However, already at the 
level of 7T° — rj mixing alone, a complete calculation at next-to- leading order produces a 
destructive interference in Qib- Additional contributions to the K — > -kit amplitudes from 
strong isospin violation at next-to- leading order essentially cancel out. The final result 
^ib — (15.9 ± 4.5) • 10~ 2 is consistent within errors with the findings of Ref. jS]. Inclusion 
of electromagnetic effects slightly increases Qib and generates sizable A and / 5 / 2 , which 
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interfere destructively with f^B to produce the final result fl e s = (6.0 ± 7.7) • 10 2 . 

It turns out that A is largely dominated by electromagnetic penguin contributions. 
In those theoretical calculations of e' where electromagnetic penguin contributions are 
explicitly included one may therefore drop A to a very good approximation. Finally, 
if all electromagnetic corrections are included in theoretical calculations of ImAo/ReAo, 
lm.A2fR.eA2 and Re/^/Re/lo, £l e s is essentially given by Qib- In this case, Q e s = (16.3 ± 
4.5) • 10~ 2 is practically identical to the result based on tt° — t] mixing only [21]. 

9 Conclusions 

In most processes isospin violation induces a small effect on physical amplitudes. In K —> 
7T7T decays, however, it is amplified by the AI = 1/2 rule: isospin breaking admixtures of the 
dominant AI = 1/2 amplitudes can generate sizable corrections to AI > 1/2 amplitudes. 
Understanding isospin violation is crucial for a quantitative analysis of the AI = 1/2 rule 
itself and for a theoretical estimate of e'. 

The theoretical description of K decays involves a delicate interplay between electro- 
weak and strong interactions in the confinement regime. Chiral perturbation theory pro- 
vides a convenient framework for a systematic low-energy expansion of the relevant am- 
plitudes. In this paper we have performed the first complete analysis of isospin violation 
in K — > 7T7T decays induced by the dominant octet operators to NLO in CHPT. We have 
reported explicit expressions for loop and counterterm amplitudes, verifying cancellation 
of ultraviolet divergences at NLO. 

On the phenomenological side, the main features/results of this work are: 

1. We have included for the first time both strong and electromagnetic isospin violation 
in a joint analysis. 

2. Nonleptonic weak amplitudes in CHPT depend on a number of low-energy constants: 
we have used leading large- estimates for those constants which cannot be obtained 
by a fit to K — > titt branching ratios (i.e., all NLO couplings and the electroweak 
coupling of order e 2 Gsp )- Uncertainties within this approach arise from (i) input 
parameters in the leading 1/N C expressions as well as from (ii) potentially large 
subleading effects in 1/N C . We have discussed the impact of both (i) and (ii) on the 
relevant quantities. 

3. Using this large- iV c input, we have performed a fit to the CP-even component of the 
couplings g 8 and g 27 , both without and with inclusion of isospin breaking. We find 
that in general the inclusion of NLO effects (loops and counterterms) has a significant 
impact on the output. The main outcome of the NLO fit is that both g 8 and #27 
are only mildly affected by isospin breaking (e.g., #27 gets shifted upwards by only 
2%). While this result is fully expected for g$, in the case of #27 it arises from non- 
trivial cancellations between LO and NLO corrections. For the ratio measuring the 
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AI = 1/2 enhancement in K° decays we find ReAo/Re^ = 20.3 ± 0.5, compared to 
22.2 ± 0.1 in the isospin limit. 

4. Using as input a NLO calculation of electromagnetic corrections to tttt scattering [55] 
155] . we have used the optical theorem to study the effect of isospin breaking on the 
final-state- interaction phases [5]. According to our analysis, isospin breaking leads 
to a discrepancy between the theoretical prediction [53] of 5 — 5 2 from pion-pion 
scattering and its phenomenological extraction from K —>■ tttc (see also Refs. [EUS]). 
Before drawing a definite conclusion about the possible presence of an additional 
AI = 5/2 amplitude, more work is necessary to understand this discrepancy. 

5. We have studied the effect of isospin violation on the direct CP violation observable 
e'. In this case isospin breaking affects the destructive interference between the two 
main contributions to e' from normal and electromagnetic penguin operators. Apart 
from the traditional term Qib, we have identified and studied the effect of isospin 
violation in the ratio ImAo/ReAo, parametrized by the quantity Ao and the purely 
electromagnetic A J = 5/2 amplitude. Both Ao and the AI = 5/2 contribution / 5 / 2 
interfere destructively with f2 IB to yield a final value fi eff = (6.0 ± 7.7) • 10~ 2 for the 
overall measure of isospin violation in e'. If electromagnetic penguin contributions 
are included in theoretical calculations of Iim4 /Reylo, A can be dropped in fl eS to 
a very good approximation. Finally, if all electromagnetic corrections are included in 
ImAo/KeAo, ImA2/ReA2 and Re/b/ReAo, f2 e ff is essentially determined by ^ib and 
is practically identical to the result based on tt° — r\ mixing only. 
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A NLO effective Lagrangians 

In this appendix we collect the relevant parts of the NLO Lagrangians. 

First we recall our notation. The gauge-covariant derivative of the matrix field U is 
denoted D^U, the external scalar field x accounts for explicit symmetry breaking through 
the quark masses, the matrix A = (A6 — iX 7 )/2 projects onto the s —> d transition and 
Q = diag(2/3, —1/3,-1/3) is the quark charge matrix. For compactness of notation, we 
introduce the definitions 

Qu = WQU. (A.l) 

Starting with the strong Lagrangian (|2.5|h we have the familiar terms JT] 

Y^UOf = L 4 (D^D»U)(x u + ) + L 5 (D^D»U X U + ) 

+ L 1 { X l i) 2 + L,{x ] Ux ] U + xU ] xU ] ) + ... (A.2) 

For the explicit form of the strong Lagrangian of 0(p 6 ) we refer to Ref. [TB*] . 
The electromagnetic Lagrangian ()2.8|) is explicitly given by [F3\ 

E K i°f p2 = K 1 {D^D^U){Q 2 ) + K 2 {D^D^U){QQ U ) 

+ K 3 ({DpUtQU) 2 + (DrUQrf) 2 )+K4(DrUtQU)(DJJQUt) 

+ K 5 {{D^U\ D^U}Q 2 } + K 6 (D^D^UQU^QU + D^UD^QUQU^) 

+ K 7 (x u + )(Q 2 ) + K 8 (x u + )(QQu) 

+ K 9 ((x^u + u^x)Q 2 + (xU 1 + u x 1 )Q 2 ) 
+ K 10 (( x j u + tf x )QtfQU + ( x u j + U X ] )QUQU^) 
+ K n (( X j U - U^ X )QU j QU + ( X U j - U X ] )QUQU^) 
+ K 12 (D^[D»Q R , Q]U + D^U[D^Q L , Q]W) 

+ K 13 (D^QrUD^QlU^) , (A.3) 
with (1^ and r M denote external spin-1 fields) 

D,Q L = d,Q-t [l„ Q] , D^Q R = d,Q-i [r„ Q] . (A.4) 

Turning to the nonleptonic weak Lagrangian, we first display the octet couplings in the 
notation of Ref. ^B] : 

E N i°i = N 5 (\{ X U + ,D^D»U})+N 6 (\D^U)(tfD»U X U + ) 

+ N 7 (X X U + ) (D^DfU) + N s (XD^D^U) ( X u + ) 

+ N 9 (X[ X U _, D^D'U]) + N 10 (A(^) 2 > + N u (X X u + )(x U + ) 

+ iV 12 (A(^) 2 )+iV 1 3(A^)(x^) + ... (A.5) 
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The corresponding 27-plet couplings [12] are 



£ AOf = £>i fii,fci (Ai iX ?) (A fc «x?> + ^2 t ijM {\ jX u _) (X klX -) 

- D i t ij , kl ()* j lflD ll U)(\ kl {x l l,rfD' M U}) 

+ D g U iM (\ ijX U + )(\mD^D»U) 

- D 7 Um (XijU^D^U) (\ kl U^D»U) (x?> (A.6) 
with (Xij) a b = S ia Sjb. The non-zero components of t^i are given by k, I — 1, 2, 3) 

^21,13 = ^13,21 = ^23,11 
^22,23 — ^23,22 = ^23,33 

Finally, the relevant part of the electroweak Lagrangian of 0(e 2 G$p 2 ) [22] is 

Y,ZiOf w = Z l {\{Q Ul x U + }) + Z 2 {\Q u )(x U + ) + Z,{\Q u ){x U + Qu) 

+ Z A {Xxl){QQu) + Z b {\D^D»U) + Z 6 (X{Qu,D^D^U}) 

+ Z 7 {\D P JJ ] D'"U){QQ U ) + Z % {\D^U){QU^D^U) (A.8) 

+ Z 9 (\D^U)(QuU^D^U) + Z 1(i {\D^U){{Q,Qu}U ] D^U) 

+ Z n (X{Qu, U ] D^U})(QD^U) + Z l2 (X{Qu, U ] D^U}){Q V D^U) + ... 

B Explicit form of NLO loop amplitudes 

In this appendix we report explicit expressions for the NLO loop corrections A L A^ ap- 
pearing in the master formulas of Eqs. ()4.6|) and (|4.7|) . 



t 



11,23 



t 



33,23 



(A.7) 



B.l Photonic amplitudes 

Let us start with the terms arising from exchange of virtual photons. The amplitude A+- is 
infrared divergent and is regulated by introducing a fictitious photon mass M 7 . Moreover, it 
is convenient to work with a subtracted amplitude, after removing the infrared component 
A 1 ^ (see discussion in Sec. HJ). The function £> + _(M 7 ) appearing in our definition of the 
infrared- divergent amplitude A l +_ (see Eq. (|4.5j) ) is given by 



B+_(M 7 



1 

47T 



M 2 1 + 3 2 1 + 3 

2a(/3) log ^ + -±jj-h{ff) + 2 + 3 log ^ 



+ ITT 



'l + (3 2 



log 



M 2 K 3 2 
M 2 
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where 



a(/3) 
/(*) 



(1 - AMl/M 2 K ) l/2 

l+£ 1-/3 
1 H log 

2/3 s 1 + /3 



1-/3 4/3 2 



2/3 



2/3 



- —log II — t| 
Jo t 



The amplitudes A^.44 7 ) are given by 

V2 



A L A\ 



(7) 

/2 



(4ttF,)2 

1 

(4^F 



H Mw 2 + 2Ml(l + log^ 



+ 



4v/2M 



*A(i/ x ) 



+ 



A L „4 



(7) 
5/2 



6 Ml-M^ / Jg\ 12(Mj-M 2 ) 

5 (4^)2 l i + l0g ^)+ 5F 2 



The divergent factor A(z/ X ) is defined in Eq. (|2.11j) . 



(B.2) 



(B-3) 



(B.4) 
(B.5) 



B.2 Non-photonic amplitudes 

The mesonic loop corrections can be expressed in terms of the following basic function 
(and its derivatives): 



J(p 2 ,Mf,M 2 ) = 

The subtraction term is given by 
J(0,M 2 ,M 2 ) = 
T(M 2 ) = 



d d k 



1 



i J (2Tt) d [k 2 - M 2 ] [{k - p) 2 - M 2 } 
J(p 2 ,MlM 2 ) + J(0,M 2 ,M 2 ) . 



M 2 T{M 2 )- M 2 T{M 2 \ 

M 2 - Ml 

1 , M 2 
log— 5- • 



" 2AK 



(4tt) 



(B.6) 

(B.7) 
(B.8) 



Expansion around the neutral meson masses generates terms involving derivatives of the 
function J(p 2 , M 2 , Mf ). In order to deal with such terms we use the notation 



J(lA0)( p 2 )M 2 )M; 



_d_ 

dp 2 



J(p 2 ,M 2 ,M 2 ) 
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J(0,l,0)(p2 >M 2 }M | 



J^ 1 \p\MlM 2 2 ) 



d 



J(j?,M*,Mi) 
J{ P \MlMl) 



(B.9) 



We report below the explicit form of the relevant functions. For this purpose we define 



\{t,x,y) 



t-{yx + ^yj t — U/x 



and 
Then 



£ 12 = M 2 + M\ 



J{p\MlMl 



32tt 2 



2 + A i 2 M l 



V 



2 log M 



A12 = M 2 - M| 

^12 



(B.10) 
(B.11) 



A 



12 



Ml 
Mf 



log a 



J{p 2 ,M 2 ,M 2 



^j^MlMl 
V 2 



2 — a log 



16tt 2 



x log 

1 



[p 2 + A 1/2 (p 2 , M 2 , M|)] 2 - A 



2 N 
12 



[p2_ A l/2 (p2)M 2 ;M 2 )]2 _ A 2^ 



(T 



^J\(l,M 2 /p 2 ,M 2 /p 2 ) . (B.12) 



The relevant derivative functions are reported below (recalling the symmetry property 

j(0,l,0)( p2)M 2 )M 2) = J(0Al)( p 2 )M | )M 2^. 



J(lA0) (p2)M 2 )M2) 



32vr 2 



A12 . Mf 

(p2)2 1Q S M 2 



(p 2 E 



12 



A? 2 ) 



(p 2 ) 2 X 1 / 2 (p 2 ,MlMl 



log 



'S 12 -p 2 -A 1 / 2 (p 2 ,M 2 ,M 2 )' 



JW)(p2,M?,A^) 



A 2 

^12 



v E 12 -p 2 + A 1 /2( p 2,M 1 2 ,M|) 
2M 2 p 2 



32tt 2 



+ 



I A 

p 2 + A 



12 



p 2 A 2 



log 



12 



m[ 

M 2 



1.2 



p2 A l/2( p2)M 2 )M |- 



log 



J(0,0,l) (p 2 5M 2 ;M 2 



32tt 2 



1 

M 2 



p 2 a 



log 



/ S 12 -p 2 -Ai/y,M 2 ,M 2 ) 
{z 12 -p 2 + \V2(p 2 ,MlMl) 

a + V 

a — 



(B.13) 



We recall here that we expand all our amplitudes around the neutral pion and kaon 
masses M n and Mk to define the isospin limit (see Ref. [BO] for a more general discussion 
of the splitting between strong and electromagnetic contributions). This applies also to the 
f] mass given in Eq. (|3.4j) . Therefore, in all (loop) amplitudes where M 2 appears explicitly 
it actually stands for (4M^ — M 2 )/3 instead of the physical value in (j3.4j) . This concerns 
all loop functions in Apps. B and C. 



In the equal mass case we adopt the definition 



j(°Al)(p2 )M 2 )M2) 



lim J^°^(p 2 ,MlM 2 

M 2 ->M 2 



1 d 



2 dM 2 



J{p 2 ,M 2 ,M 2 ) . 



45 



B.2.1 AI = 1/2 amplitudes 

In this section we list the one-loop corrections to the AI = 1/2 amplitude. 



A,4/? = -g| J(M|, Ml Ml) + (2 Ml) J(M|, Ml Ml) 

+ Ml Ml - M « ( gj ^ M ^ J ( M 2 , M|, M 2 ) 

3(12Mj-llM|M, 2 + 3M, 4 ) J 2 * (6M|-11M|M 2 ) 

8F2(-M| + M2) 1 4F2(M1-M2) 1 ^ J 

(8 M| — 35 M\ Ml + 25 M*) 2 -M| + Ml 
+ 8 Ft (M* K -Ml) T(MJ+ 16 F*** (R14) 

= M, 2 , M;) + (2 ^ " J(M|, Ml Ml) 

' f ' [ AMI Ml, Ml - {M4k ;^ M J MI) J {Ml Ml, Ml) 



YlFZMl v " y 1 /•- Ml 



2 



(36 M| - 73 M| M 2 + 19 Mp M K 

72 F* (Ml -Ml) 1[M v ] + AFl 1 {Mk} 

(8MI-35MIMI + 25MD -9A# + 4Aff 

+ SFjjMl^Mj) i(MJ+ 144 Ff^ (B - 15) 

+ (5 M \ '.I K) AMI, Ml Ml + gjg^jg) J(A4, M-) 

- (3 M| -4 M| M?) 2 M2 2 _ (Mj - 2 M| M 2 ) 2 2 ^ 
\2F 2 M 2 v F 2 M 2 

(22 M| -71 M| M 2 + 43 M*) 5M|M 2 

12 F2 (Ml -Ml) 1 " ;_h 4F2 (Ml -Ml) 1 ^ J 

(4M|-11M1M, 2 ) 8M ; - 22 M^M: • 7 M: 1 



4F 2 (M1-M 2 ) 32 F 2 (4 Ml -M 2 )^ 

r2 /■ j\/f2 



+ 



4M-K-M-) j(0 , , 1)(m , m , m;) 



9^ 2 
Ml (Ml - M 2 ) 



3 F 2 Ml 

7T 7T 



J^\Ml,Ml,Ml + /^(M^Ml^ 2 )] (B.16) 



AlA[ Z /1 = -^t «, M|, M|) - (2 M | ^ M - } J(M|, Ml Ml) 



(2MI + 5MIMI-4MIMI) 

2AF 2 Mt J(M W ,M K ,M V ) 



TT TT 
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(Ml -A M 2 K Mi) -, * o 9n 

1 \ F 2 M , ^ mi mi, md 

Ml (A Ml -Mp _ 3(7Ml-7MlMl + Mt) 2 

8F2(Ml-M%) 1 " J 8F2(Mi-M2) 1 ^ J 

(7M|-18M2) 8j4-13M|M^+2^ 

8F2 ^ 128F2M2 7T 2 

- < 2 ^ - 8 f 2 M ' + J^(Mh Ml Ml) 



+ 
+ 



(Ml - 5Ml Ml + AMI M$) 
A Ml 



J(0,0,1) (M 2 ? M 2 } + Jd,0,0) (M 2 ? M 2 ^ M 2^ 



(M| M 4 K Ml) rj(o,i,o) (M 2 ;M | )M 2 ) + 7(1.0,0)^2^^)1 (B17) 



12F2M2 
r2 



= Mi) + (2 M ^ 2 MJ J« M2, M^) 



7T 



(8 III - 6 Ml Ml + Ml) | ^Mj^MjM, 2 ) 

8F2(M|-M2) 1 + 8F2(Mi-M2) 1 ^ J 

(8 Mi - 35 M 2 , Af2 + 21 AT*) , 9 , -5MU4M^ 
+ 8^(4-^) T(M ' )+ 128V • <B ' 18) 

B.2.2 A/ = 3/2 amplitudes 

In this section we list the one-loop corrections to the AI = 3/2 amplitude. 



= - (M| 2 / 2 M ' 2) ^ *0 ~ ^^2 J (Ml Ml, Ml) 

TV TV TV 

Ml (5 Ml -8 Ml) ... - „, 

~8fJmi — J(M '' M ' } 

M2(4M2- ?r M2 ) (3 Ml + Ml Ml) 

+ 8F2 (M| - M 2 ) Km *> + 4F2 (Mi - Ml) { K> 

(AMI -22MI Ml + 29 Ml) , Mj-2Mj 
8^ 2 (M|-M2) T(MJ+ 32F2.2 (B.19) 

A*4$ = +7^| ^, M*) - (Mi ^ M - } J(M|, M 2 , M 2 ) 

_ (21Mi-8MjM2) * _ (llMj-16MjM|) 

24F 2 M 2 J \ M ^ M K, M r,) 8F 2 M 2 J V M n, M n, M K) 



TV TV TV 
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(44 M| - 47 Ml Ml + 9 M 4 ) (12 M 4 - 17 Mj M 2 ) 

24F 2 (M£-M 2 ) 1 " J 4F 2 (M|-M 2 ) 1 k) 

(AMI - 11 Ml Ml + 15 M^) | 2M|--M 2 

8F2(M|-M2) 1 ^ 16F 2 7r 2 

+ ^^^ J( °' 1 '°H^^^^ 2 ) (B-20) 

(10 M| + 13M| M 2 - 32 Ml M 4 + 24 M®) 2 2 2 
12M 



7T 7T 



(lOMj + 3M 4 M 2 -28MlM 4 ) 

40F2M 4 -W, MJ 

(48 M 4 - 40 Ml Ml + 7 M 4 ) 3 (21 M 4 - 20 M| M 2 ) 

40F 2 (M£-M 2 ) 1 " J 20F 2 (M£-M 2 ) 1 k) 

(58 M 4 - 22 Ml Ml - 27 M 4 ) 2 -M 4 + 14 M| M 2 - 10 M 4 

40 F 2 (Mi - Ml) [ 80 F 2 M 2 tt 2 

2 (M£ — 3 Ml Ml + 2 M 4 ) moi)/lijr2 A/r2 A/r2 , 

+ -A— - — 5 ^ 2 ^/"-^(A^XX) 



+ 



+ 



(5M| - 13M 4 Ml + 8M 2 M 4 ) 

20 F 2 M 2 
(M|-M 4 M 2 ) 



J(1,0,0) (M 2 ; M 2) + J(0,0,D (M 2 ) M 2_ ; M 2) 



12 F 2 M 2 



J(°,m(Ml,Ml,M*) + /^(M^Ml^^ 2 ) (B.21) 



7T 

_ (5M 4 -8MjM 2 ) * 

8F 2 M 2 J y 1V1 n, 1V1 K, M 7r) 

(8 Ml -6 Ml Ml + M*) (2MI-5MIMI) 

8FKMI-MD 1 " j 4F 2 (M|-M 2 ) 1 * j 
(4M 4 +2M 2 -M 2 -3M 4 ) M| - 2 M 2 
8F. 2 (Mi-M 2 ) T(MJ + 32F 2 vr 2 • (B ' 22) 

B.2.3 AI = 5/2 amplitudes 

In this section we report the one-loop AI = 5/2 amplitude generated by insertions of e 2 p° 
vertices from C c \ m . 
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2 (M A K + M 2 K Ml-2M*) 

15F 2 M 2 
2 {AM 2 K -Ml) 



J(M 2 , M 2 K , M 2 ) 



5F 2 



T(M 2 ) - 



AM 4 K 



"' 5F 2 (M 2 K -M 2 ) 
T(M 2 ) + - 



T{M 2 K ) 



2 (6 Mr — 19 M 2 K Ml + 11 Ml) 2 -Af£ + 9 M£ - 10 M* 



5F2(M|-M2) 

4 (Mi-3M^M2 + 2M4) - (001) . „ n/r2N 
+ ^J mi) (Ml,M 2 : M 2 ) 

7T 

4 (Mi - Ml M^ 



40 F2 M2 7T 2 



5F2 



B.2.4 Divergent parts 

For completeness, we list here the divergent parts of the mesonic loop amplitudes. We have 
checked explicitly that they get absorbed by the independently known renormalization of 
NLO chiral couplings. 



(e) 
1/2 

(zy 

1/2. 



A L A 



A L A 



\ L Af> 



(e) 
3/2 



A L A 



A L A 



(ZY 
3/2 



A 



ArA {Z) ' 



div 



div 



div 



div 



div 



div 



div 



div 



div 



div 



-28M|. + 17M2 
2F 2 



-27MI + 103M2 
18F 2 



10MI-43M2 
7(M 2 K + M 2 ) 



2F 2 
-M 2 K + 10 M 2 
2~F^ 

7T 

-MI-15M2 

2F2 
64M|-27M2 

17 (4 M 2 K + M 2 [ 



10 F 2 

SM^ + M 2 



A("x) 
*A(i/ x ) 

A(^x) 
AK) 



2F2 
48 (M£ - M2) 



5F2 



A(^x) 
A(^x) • 



(B.24) 
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C Alternative convention for LO LECs 



In the effective chiral Lagrangians of Sec. 12.21 the meson decay constant in the chiral limit 
F is the only dimensionful parameter in addition to the Fermi coupling constant Gp. This 
is the original convention of Cronin [14 for the nonleptonic weak Lagrangian of lowest order 
and it is used throughout this work. It has definite advantages, e.g., for the renormalization 
of the various Lagrangians. 

However, this convention has a certain aesthetic drawback in that the K ^ 2n am- 
plitudes (the K — > 3tt amplitudes as well, for that matter) depend at NLO on the strong 
LECs L4 and L5 even in the isospin limit. These LECs account for the renormalization of 
F to F n and F K at NLO. The associated uncertainties propagate into the uncertainties of 
the LO LECs G$, . . . Since F n and Fjc are much better known than F, it may be useful 
for phenomenological purposes to redefine the LO LECs so that they are then free of the 
uncertainties in L\,L\. 

A first step consists in generalizing the convention first used in Ref. [15] . albeit with a 
different notation: 



Gs — G 8 F 4 / F* , 

9ewk dewk-F / F^ , 



G27 — GrjF A I F 4 
2 = ZF 2 /F 2 . 



(CI) 



At lowest order, the barred quantities are identical to the original unbarred ones because 
we always set F = F n at lowest order. Writing the NLO amplitudes ()4.6|) in terms of the 
barred LECs of lowest order, the strong LEC L\ disappears completely from all K — > 2tt 
amplitudes. To get rid of L5 as well (at least in the isospin limit), one can introduce a 
scale factor F^/Fk [SHI- The amplitudes of Eq. (j4.fi|) then take the following form: 



where 



An = g 27 fJm 2 k ~mAa^ 



(C.2) 



+ G 8 F n 




A^ + e&> A^ 



e 2 F 2 



A^ ~\~ Z A^ n ^ + (^ewk A^ 1 



A^ 



Ft, 



a™ 

n T? 



K 



1 + A L AY ] + A c A { n x) ] if af ) ^ 



(C.3) 



A L AiP + A c A^ if af)=0. 

The change in notation only affects those amplitudes that are non-zero at lowest order. 

The amplitudes A^ are related to the original A^ as follows (only amplitudes for 
n = 1/2 or 3/2 are affected): 



for X = 27, 8, e : 



AcA^ = A c A^\ c=Lr _,+ 



24(M| - Ml) 
F 2 



L r 5 (v x )5 n< i/ 2 5x, l 
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A L A^ = A L ^ + A K + 3A,-^(E K + 3E^)5 ntl/2 S X!E] (C.4) 



for X = Z, g : 



AcA^ = Acv4W|~ r)=0 

A L AW = A l A^ + ^k + 5A w . (C.5) 



From the definitions of F n and Fk± in Ref. [23 one obtains [T(M 2 ) is defined in flB.8 
and M 2 stands for (4M| — M 2 )/3 as in all loop amplitudes] 

A. = ^T{Ml) + ^T{Ml) 



K = - y -^—^T{M 2 K ) -\- 



E 



K 



7T 7T TT 

+ (4^ • (C - 6) 

As can be seen from Eqs. (10.21 . . . JC.6|) . L4 has disappeared completely from the ampli- 
tudes ,A n whereas L5 occurs only in the isospin violating amplitude A c Ay 2 . In spite of its 
conceptual advantages, we have not used this alternative convention in this paper because 
L 4 , L 5 reappear anyway through the large- N c relations for the NLO LECs iVj, Di, Zi. More- 
over, the large- N c relations for g 8 , g 2 7 and g e wk would also be affected. Finally, consistent 
with the expansion to leading order in 1/N C , L\ and L r 5 are set equal to their large- N c 
limits as discussed in Sec. 15.31 



D Details on the optical theorem analysis 

In this appendix we report the explicit form of functions needed when studying the unitarity 
condition in the presence of isospin breaking. Let us start with the IR divergent factors: 



B + _(M, 
16vr 



(see Eq.flDJ) 
-2M 2 )G+_>) " (t 



2Ml) G+_ 7 (t) 



(D.l) 
(D.2) 



The definition of the variables t, u and the function G + _ 7 (x) can be found in Refs. [5*5^ 156) . 

In order to define the remaining ingredients, we need to fix the notation. The four- 
momenta are denoted as follows: 



K (P) ^ vr +(g + ) vr -(g_) 7 (A;) ^ vr +( p+ ) vr -(p_) . 
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The differential phase space is then given by 

= wm wm> {2 * } ' * (4)(?+ - . < M > 

Then, after performing the sum over photon polarizations, the radiative amplitudes in 
leading Low approximation generate the following factors: 

(D.4) 



/■rad 
Jl 








+ 


2g+ • q- 


(?+• 








rva,d 
J 2 
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p- ■ q- 


ip+ ■ k 


- -f)m 


+ ' (p_ 




_Jfi )( ,_. k + % 
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(p+ ■ k 




•* + ¥) (p- 


• fc 


M l\l 1 I M l\ ' 



(D.5) 
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